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1. Introduction 



In this paper, we consider the XX-chain with diagonal and non-diagonal boundary terms: 
1 L_1 1 

H = - J2 [vfvj+i + vjvf+i] + [a-^r + a + a t + «W + #f ff L + P-*Z + /w] ( L1 ) 

Here, are defined by o" 1 * 1 = \{p x ± id 5 '), where a x , a v and <r 2 are the Pauli matrices. The factor 
has been introduced for later convenience. Since the parameters a± , (3± , a z and f3 z are arbitrary 
complex numbers, the Hamiltonian defined by eq. Ql,l| ) is non-hermitian in the general case. 

Let us now give a brief overview of the literature before turning to the concrete results we 
obtained by studying diagonal and non-diagonal boundary conditions. All the articles mentioned in 
this overview are based on the free fermion approach to the XX-model. 

The XX-model often appears as a special case of the XY-model. The XY-model was introduced 
1961 by Lieb, Schultz and Mattis [[jj who computed its ground state energy, the elementary 
excitations and also presented a method to calculate time—independent correlation functions. In 
this way, they treated periodic boundary conditions as well as free ends. 

During the last thirty years, the correlation functions of the XY-model and therewith the XX- 
model have been the subject of various investigations. McCoy |2| studied the correlation functions 
of the XY-model with periodic boundary conditions. More precisely, he computed the asymptotic 
behaviour of each of the three time-independent correlation functions (oo "}?) with i = x, y, z in the 
limit R — > co. Barouch and McCoy || determined the asymptotic behaviour of the same correlation 
functions for the XY-model with an external time-independent magnetic field in the z-direction. 
In another article, time-dependent spin-spin correlation functions of the form {<j x (i)<r^(0)) and 
(<7q (t)o"p(0)) for the XY-model in an external magnetic field again in the z-direction were calculated 
in the limit of large R by McCoy, Barouch and Abraham ||] . Exact expressions for these correlation 
functions for all values of R and t were then computed by Vaida and Tracy || • Furthermore, time- 
dependent many-spin correlation functions for the XY-model in an external constant magnetic field 
in the z-direction were treated by Bariev k|. 

Recently, the XY-model with boundary terms has been a subject of rising interest. Hinrichsen 
and Rittenberg |7| showed that the anisotropic XY-model in an external magnetic field with o z - 
boundary terms is invariant under certain quantum group transformations. Furthermore, they 
defined and calculated the corresponding invariant correlation functions. 

The XX-model with non-diagonal boundary terms, however, has not been studied thoroughly 
before. Some work in this direction has been presented by Guinea Q who studied the semi-infinite 
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XY-model with one cr^-boundary term (i. e. a_ = q + = 1, (3 + = /?_ = a z = (3 Z = 0). We will 
mention more details of that paper when discussing the physical applications of the Hamiltonian H 



of eq. Furthermore, a study of the totally asymmetric AX-model with bulk terms of the form 

<Tj~aJ +l and with boundary parameters given by a_ 7^ 0, (3+ 7^ 0, a + = /3_ = a z = j3 z = in the 
notation of eq. ( |1.1| ) can be found in our previous paper 

As already mentioned in the beginning, in the general case, the Hamiltonian given by eq. (|0|) is 
non-hermitian. Interesting physical problems involving non-hermitian Hamiltonians can be found 
in several articles treating non-hermitian quantum mechanics |10|. 



The Hamiltonian given by eq. ( |1.1[ ) can also be used in the study of asymmetric bulk terms. 
More precisely, starting from a Hamiltonian of the form 

L-l 



H 



+ 



1 
71 



a'_ a l + a' + + a z a{ + (3+ a\ + f3'_a L + j3 z a z L 



(1.2) 



one can use a similarity transformation (see for example jjllj] ) that transforms the asymmetric bulk 
terms depending on the two parameters p and q of the Hamiltonian given by eq. ( |l-2j ) into symmetric 
bulk terms. It is convenient to choose ^/pq = ^. Note that the similarity transformation changes 
the boundary terms. The corresponding transformed boundary parameters a_,a+,/3+ and [3- of 
the Hamiltonian given by eq. (1.1) are L-dependent now and have the expressions 



Q—a'_ 



a a 



(1.3) 



with Q = The diagonal boundary terms remain unchanged. 

Although the Hamiltonian of the XYZ-chain with non-diagonal boundary terms is known to 
be integrable [12, 13], Bethe-Ansatz equations have not yet been obtained, because it is not clear 
how to construct a reference state. Therefore, to study the effect of non-diagonal boundary terms, 



we chose the XX-model with boundaries of the form given by eq. (bl), because this model can be 
fermionized. To be able to use the free fermion approach, we introduce a new Hamiltonian Hi ong 
which is bilinear in terms of fermionic creation and annihilation operators. This approach has the 
major advantage that we have complete control over the wave functions for a large class of boundary 
parameters which enables us to calculate correlation functions. Thus, we get a good handle on a 
particularly interesting and important integrable model. 

As mentioned above, in order to treat the Hamiltonian given by eq. ( |1.1[) we transfer the 
diagonalization problem to a new Hamiltonian which we obtain by appending one additional site at 



each end of the chain as in [14]. This new Hamiltonian has the following expression 



L-l 1 



H long = J E H^j+l + °j 4+1 
J'=l 



^8 

+ f3 + a+al +1 + (3. ala x L+1 + (3 Z a*\ (1.4) 



In this way, the boundary terms are also bilinear expressions in the a + and a~ matrices. It is 
only after this transformation that we can write and solve the problem in terms of free fermions. 
Since erg and <r£ +1 commute with Hi ong , the spectrum of Hi ong decomposes into four sectors 

(++, H — , — h, ) corresponding to the eigenvalues ±1 of erg and <r£ +1 . The original Hamiltonian 

corresponds to the (++)-sector. A substantial part of this paper is devoted to showing how the 
eigenvalues of H are obtained by projecting onto this sector. 

The Hamiltonian Hi ong which we introduced only as a means to treat the Hamiltonian H is 
actually interesting in its own right as a quantum spin chain with boundary terms. 

In the field-theoretic approach the Hamiltonian Hi ong is probably related to the decoupling point 
of the boundary sine-Gordon model. The corresponding boundary S-matrix has been calculated in 



0,0. 

It is very likely that the Hamiltonian H given by eq. ( |1.1[ ) can be applied to physical problems, 
since a simpler version of this Hamiltonian has already found such applications. Namely, the semi- 
infinite XX-chain with one cr^-boundary term mentioned before was studied by Guinea || as a 
model for the dynamics of a particle in an external potential coupled to a dissipative environment. 
He also utilized free fermions and presented an explicit solution for the mobility of the particle in the 
continuum limit. Afterwards this solution was used in the study of transmission through resonant 
barriers and resonant tunnelling in an interacting one-dimensional electron gas, cf. Kane and Fisher 



[ IT | . This type of system is studied in experiments with quantum wires. The calculation is built on a 
perturbative renormalization group analysis in different limits (limits of a weak barrier and a strong 
barrier). By combining the results of these two limits the authors obtain the full phase diagram of 
the model. For one particular value of the dimensionless conductance, they even obtain an exact 
solution for the conductance through a resonance by mapping the model onto the semi-infinite 
XX-model with one <7 x -boundary term. 

The starting point for our investigations of the XX-chain with non-diagonal boundary terms 



is the diagonalization of the Hamiltonian H (eq. (1.1)). This problem is not only of mathematical 
interest, since the model has an interesting physical content. Namely, as will be shown, boundary 
bound states appear and the non-trivial ground state expectation values of the aj-operators and 
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the crj-operators exhibit a decay into the bulk which can be predicted from conformal field theory. 
Furthermore, the expressions for the partition functions formally coincide with partition functions 
of a Coulomb gas with only magnetic charges or only electric charges, depending on the choice of 
the boundary parameters in the Hamiltonian H. Additionally, the fermionic energies as well as the 
expressions for the ground state energies show a logarithmic dependence on the lattice length for 
special choices of the boundary parameters. The study of these physical properties is deferred to two 
subsequent articles. Below we summarize the content of all articles and point out how the results of 
the present article enter into the further considerations. 

In this first article, we confine ourselves to studying the integrable model with non-diagonal 
boundary terms given by H on a finite chain. This includes the calculation of the spectrum and 
the wave functions as well as the derivation of expressions for the one- and two-point correlation 
functions for the crj-operators where the subscript j indicates the position on the chain. These 
results are obtained in parallel for H and Hi ong . Let us briefly describe how we proceed: We start 
by fermionizing the Hamiltonian Hi ong . The spectrum of the original chain as well as the eigenvectors 
can be retrieved from the spectrum and the eigenvectors of Hi ong by a projection technique which we 
derive in detail. As a by-product, we solve the eigenvalue problem for the quantum spin chain H.\ ong . 
We demonstrate that after fermionization the problem of finding the eigenvalues of the Hamiltonian 
Hiong is reduced to the problem of finding the zeros of a complex polynomial of degree 2L + 4, which 
we write down explicitly. This polynomial, which might very well also appear in other contexts, has 
interesting algebraical properties. Namely, for some choices of the parameters, it can be factorized 
into cyclotomic polynomials. We looked systematically for these factorizations since, apart from 
being of mathematical interest, these examples give access to an exact solution for the full spectrum 
of the Hamiltonian, including exact expressions for the ground state energy. (In the general, non- 
hermitian case, we define the ground state energy to be the one with the smallest real part.) Some 
of these examples are especially interesting since the factorizations contain L-independent factors 
which lead to L-independent eigenvalues of the Hamiltonian. The corresponding eigenstates will 
be identified as boundary bound states in the next paper. Furthermore, the "cyclotomic" examples 
furnish a reliable ansatz for an approximative study of the zeros of the polynomial in the general 
case which will be presented in the third paper. As an additional result, we get exact formulas for 
the one- and two-point correlation functions for the crj-operators. The value of enters the 

projection mechanism mentioned above. 
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In the second article, by using the results of the first paper, we calculate one-point functions 
for the a?— and the crj-operators for arbitrary position j and lattice length L for several of the 
"exactly solvable" cases where the polynomial can be factorized into cyclotomic polynomials. These 
one-point functions decay into the bulk with a power law typical of conformally invariant theories. 
Taking this point of view, we determine their critical exponents. 

Furthermore, we make the connection between excitations with an L-independent energy seen 
in this paper and boundary bound states. This identification is made on the one hand by studying 
the spatial profile of the special fermionic excitations in comparison to the spatial profile of other 
fermionic excitations, and on the other hand by comparing them to boundary bound states found 



in the Bethe ansatz for the XXZ-chain with diagonal boundary terms [19|. Boundary bound states 
originally have appeared in the field-theoretic approach to the sine-Gordon model with boundary 
interaction jl8j, [l^, In our case, it is surprising that they are related to special zeros of the 
complex polynomial as mentioned above and can therefore be found without invoking the field 
theory. 

A further important new observation is related to the partition functions in the thermodynamic 
limit. They will be presented in the third article, where they will be derived by studying 
approximative solutions of the polynomial equation for large values of L as mentioned above. The 
partition functions correspond to conformally invariant systems, a behaviour which we also found in 
our previous study of the totally asymmetric XX-chain with non-diagonal boundary terms ||. This 
observation is confirmed by the expansions of the exactly calculated ground state energies for large 
L. From this expansion one can read off the conformal charge c = 1 and obtain expressions for the 
surface free energy. Moreover, the partition functions we find are the partition functions of a Coulomb 
gas with only magnetic charges or only electric charges. The phenomenon of finding only magnetic 
charges is elucidated by the construction of a pseudoscalar magnetic charge operator from the 
fermionic number operators which commutes with the Hamiltonian for finite chains. Furthermore, 
for special choices of the boundary parameters, we find a logarithmic L-dependence for the fermionic 
energies as well as for the expression for the ground state energies. This may only happen if the 
Hamiltonian is non-hermitian. 



The present article is very technical by nature. For those readers who are not interested in all 
details of our calculations but would nevertheless like to use our results without reading the whole 
paper we provide a guide in section 13, which does not, however, follow the sections in a chronological 
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way. The others sections are organized as follows: In section 2, we use the fermionization of the 
chain Hi ong to reduce the eigenvalue problem of this Hamiltonian to the eigenvalue problem of a 
matrix M of dimension (2L + 4) x (2L + 4) whose eigenvalues correspond to the fermionic energies. 
We derive some general properties of the eigenvectors of M (which will be needed in the sections 10 
and 11) before showing, in section 3, how the eigenvectors corresponding to the non-zero eigenvalues 
of M can be calculated explicitly. The solution of the eigenvalue problem of M leads to a complex 
polynomial (which corresponds to the characteristic polynomial of M) whose zeros determine all 
eigenvalues and eigenvectors of M. This polynomial is presented in section 4. Section 5 is devoted 
to the study of the factorization properties of this polynomial. By constraining the total number of 
cyclotomic factors, we systematically determined the boundary parameters for which the polynomial 
factorizes into cyclotomic polynomials. Some of these cases are actually one-parameter families of 
solutions. In section 6, we show for two examples how the full spectrum of Hi ong is obtained from 
the factorized form of the polynomial. Section 7 contains the exact expressions for the ground state 
energies of all examples where the polynomial factorizes into cyclotomic polynomials. In section 
8, we present one example of a Hamiltonian with asymmetric bulk terms where it is also possible 
to calculate the full spectrum and the ground state energy exactly for arbitrary values of L. In 
section 9, we derive the projection mechanism which is needed to obtain the spectrum of the original 
Hamiltonian from Hi ong . To derive the projection mechanism we need the value of the one-point 
function of the uj-operator at the point j = L + 1. We express the one- and two-point correlation 
functions of crj in terms of Pfaffians in section 10. In the cases where the Hamiltonian H has no a z 
boundary terms or fulfils the condition a_ = a + and f3 + = j3- , we further reduce these Pfaffians 
to determinants of a certain matrix. These expressions will be needed for the calculation of spatial 
profiles in our second article. In section 11, we calculate the above mentioned value of the one-point 
function of <rj at the point j = L + 1 in the cases where the Hamiltonian is: a) hermitian, b) has no 
a z boundary terms or c) fulfils the condition q_ = a + and /?+ = /?_. Inputting this result we invoke 
the projection mechanism and present the ground state energies for the original Hamiltonian H in 
the "exactly solvable" cases which additionally satisfy at least one of the afore mentioned conditions 
a)-c) in section 12. We conclude this article with a discussion of our results in section 14. In an 
appendix we show how to find the eigenvectors of the matrix M corresponding to the eigenvalue zero. 
We derive the conditions for the appearance of zero modes in the spectrum of H[ ong and determine 
respective restrictions for the boundary parameters. 



6 



2. Diagonalization of the Hamiltonian 



In this section, we present the general formalism we use for the diagonalization of the XX-model 
with boundary terms defined by eq. H can be diagonalized in terms of free fermions if it 

can be written as a bilinear expression in u^-matrices, since standard fermionization techniques can 
then be applied @3- 

To obtain a bilinear expression in cr^-matrices for H we add one lattice site at each end of the 



chain, site and site L + l as in [14|. Notice that the terms containing a z do not have to be changed. 

The Hamiltonian now reads 

L-l 



Hlong = ^ E K +(J 7+1 + a 7 a Hl 



1 



a-cr a 1 + a + a a{ + a z a l 
+ (3+ a+a% +1 + /?_ a~ L al +1 + (3 Z a z L ] (2.1) 



As Cq and cr£ +1 commute with Hi ong , the spectrum of Hi ong decomposes into four sectors 

(++, H — , — h, ) corresponding to the eigenvalues ±1 of o"g and <t£ +1 . Notice that the projection of 

Hi on g onto a fixed sector (ei, £2) has the same eigenvalues as the original Hamiltonian with the choice 
of the parameters t\ a__, ei aq_, 62 /3_ and €2 (3 + so that by diagonalizing H[ ong one simultanously 
treats four different Hamiltonians H. The eigenvectors of the original choice of the parameters 
a_,a + ,/3_ and j3 + can be retrieved by projecting onto the (++)-sector as described in section 9. 

Furthermore notice that the (++) sector and the ( ) sector respectively the (H — ) sector 

and the ( — h) sector can be interchanged by using the following transformation which leaves H[ ong 
invariant: 

< • < o) ■ a) a j -a) j = 0, . . . , L + 1 (2.2) 

It maps any eigenvector |^) of Hi ong from the (ei, 62)-sector onto an eigenvector |*)' of H\ ong with 
the same eigenvalue lying in the sector (— ei,— 62). Therefore each eigenvalue of Hi ong is at least 
twofold degenerate. In the fermionic language, the above symmetry manifests itself as a zero mode. 

In the next section, we will show that the diagonalization of Hi ong can be reduced to finding 
the eigenvalues and the eigenvectors of a (2L + 4) x (2L + 4)-matrix which will be denoted by M. 
After studying general properties of the eigenvectors, we will describe in section 3 how they can be 
obtained in an explicit form. The eigenvectors and the eigenvalues of the matrix M are determined 
by the zeros of a polynomial which will be given in section 4. 
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2.1. Diagonalization of 

Adopting the Majorana representation of the lattice s = 1/2 spin operators as in set 

rp (EK>r ( 2 - 3 ) 

i=0 

These operators obey the anticommutation relations of a Clifford-algebra {t^,t^} = 25^. 
Rewriting Hi ong in terms of Tj ~ , we obtain the following bilinear expression 

Hlon 9 = - E E^^+l + ^^ + ^^i+^^ + J^T^ (2.4) 

fM,V=±l j=l 



where 

/ i 




(a_-Qt+) 4=(a_ + a + ) \ W -j=(p+ + (3.) 




Here we chose the basis such that the matrices above have the general form A 



(2.5) 



A A-+ 

A+- A++ 

where A is one of the matrices in (2.5). Now we apply a linear transformation to the Tj ~ operators 
to obtain another set T+,T~ of Clifford operators again satisfying 

{T£,TZ} = 26% 1 (2.6) 

Let 

L+l 

n = E E (2.7) 

3=0 n=±l 

be the explicit form of this linear transformation where 7 = ±1. One can choose this linear 
transformation in such a way that in terms of these new Clifford operators Hi ong takes the simple 
form 

L+l 

Hiong = E A " iT n T " (2.8) 
n=0 

The commutation relations for the T~,T+ imply that the operator iT~T+ has eigenvalues ±1 so that 
the spectrum of H\ on g is given by all possible linear combinations involving all A n with coefficients 
+1 or —1 and can be read off eq. ( pEq) . 



Notice that the operators T n ,T+ as defined by eq. (2.7) are in general non-hermitian. However, 
according to a general theorem for Clifford operators it is possible to apply a similarity 



transformation to the set of vectors (ipn) ano - (V'n) ^° obtain new hermitian Clifford operators 
T~',T^' in terms of which the Hamiltonian also takes the form given by eq. ( |2.8| ). This will be 
discussed in detail in the next section. 

The coefficients (ip%)j of eq. are constrained by requiring that the operators obey the 
anticommutation relations of eq. ( |2.6D : By computing the commutator [Hi ong , T^] using for Hi ong 
first the expression (|2.8|) , and then ( |2.4|) , and comparing both results, one finds that the eigenvalues 
A„ and the vectors 



are given by the solutions of the following equations 
Mip+ = -LM£ , Af^~ = iA^+ 
where M is a (2L + 4) x (2L + 4) matrix given by 

G 



(2.9) 



(2.10) 



M 



-G 1 



21 F 
-F T 



—F T 



V 

Defining 
leads to the eigenvalue problem 



2J 



J 



(2.11) 



M0; 



± 



±K4>t 



(2.12) 



(2.13) 



Observe that M has 2L + 4 eigenvalues although Hi ong has only length L + 2. This can be 
explained by considering eq. ( |2.13| ). As one can see, with the appearance of each eigenvalue A n we 
also get the negative eigenvalue — A n . As mentioned above, the spectrum of Hi ong is given by all 
linear combinations of A n with coefficients ±1 (see eq. ( p.8[) ) and thus can be retrieved from the 
eigenvalues of M by choosing from each pair of eigenvalues ±A n one value as basis element for the 
Z2 - linear combinations. Later we will make this choice in a systematic way following a physical 
convention which consists of choosing as relevant energies the eigenvalues with positive real part. 

As can be seen directly from the form of M, at least two of the eigenvalues A n are zero. The 
corresponding eigenvectors are given by (0, 1, 0, 0, . . . , 0) and (0, 0, ... , 0, 1,0). Since the eigenvalues 
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of M occur in pairs ±A n , from which only one value has to be taken, the zero eigenvalues lead to one 
zero mode as already mentioned above. As we are going to see more explicitly, this zero mode does 
not appear as a fermionic excitation in the spectrum of H. This fact can be explained as follows: 
Recall that in the case of Hi ong the zero mode reflects the presence of the symmetry given by eq. 

(|2.2|) which interchanges the (++) sector and the ( ) sector respectively the (H — ) sector and the 

(— +) sector. Since H corresponds only to the (++) sector, it is clear that the above symmetry is 
not a symmetry of H. Therefore the above zero mode does not appear in the spectrum of H. In the 
following we are going to call it spurious zero mode. 

To express the spectrum of Hi ong in terms of free fermions, we will now write the expression for 
the Hamiltonian in terms of fermionic operators b n and a n satisfying 

{K, Cm} = 5 n ,m , {K, b m } = , {a n , a m } = ( 2 -14) 
which are obtained from the Clifford operators T+ and T~ by the following transformation: 

b n = \{T+ + iT-); a n = ±(T+-iT-) (2.15) 
Hi ong then reads 

L+l L+l L+l 

Hi ong = 2 A " b ™ a ™ - Yl A ™ = £ 2 A « Nn + E ° ( 2 - 16 ) 

n=0 n=0 n=0 

where Eq is the ground state energy of the system and N n the number operator (with eigenvalues 
and 1) for the fermion with energy 2A n . 

Notice that in the expression for the number operator N n in eq. ( p. 16 ) b n is equal to a n if the 
operators and T~ are hermitian. As mentioned above, they always can be chosen to be hermitian 
by applying a similarity transformation to the vectors (V'n ) an d {^n ) m e< l- ( P-7| ). At the same time, 
the operators a n and b n are then transformed into new operators a' n and b' n which are adjoints of 
each other. 

In eq. ( 2.16| ) we have defined the Fermi sea by summing over all negative eigenvalues of M. 
Consequently, we have to choose the other half of the eigenvalues of M to form fermionic excitations 
above the Fermi sea. Here and in the following we will use the convention that if a pair of eigenvalues 
has non- vanishing real part, we will denote the one with positive real part by A n . This choice leads 
to a ground state energy with the smallest real part. In the case where the real part (but not 
the imaginary part) of A n is zero one has the freedom of choice to take either the eigenvalue with 
positive or the eigenvalue with negative imaginary part as a fermionic excitation above the Fermi 
sea. This leads to an ambiguity in the value of the imaginary part of the ground state energy. A 
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similar problem occurs in the calculations involving the eigenvectors of the zero modes (e. g. in the 
calculation of one-point functions of cr-operators). Namely, the zero eigenvalues of M also occur in 
pairs C+0' and '—0') and one can freely choose which of these two zero eigenvalues belongs to the 
Fermi sea and which one corresponds to an excitation with zero energy. In other words, one can 
choose which is the creation and which the annihilation operator corresponding to the fermion with 
zero energy. We will come back to this point in [p^ ]. 

The fermionic operators a m , b m can be expressed in terms of the Tj, tJ -operators by using the 



eigenvectors of M and eq. (2.7) 



1 1 L+1 

am = 2 " = - £ £ (4> m )W (2.17) 

j=0 fj=±l 

bm = \{T+ + IT') = \Y, E ^T 3 rf (2.18) 

Notice again that in general the expression given by eq. ( p,18| ) is not the adjoint of a m because 
the Hamiltonian is non-hermitian. However, if one transforms the operators T n ,T+ into hermitian 
operators as mentioned above, the set of vectors ((fr m ), {<t> m ) fulfils the conditions (4> m ) = (<f> m )* and 
b m becomes the adjoint of a m . 

2.2. Orthogonality relations 

In the following, we make some general remarks on the given eigenvalue problem for the skew- 
symmetric matrix M = —M l defined by eq.( 2.13| ). We show that we can indeed find a linear 



transformation of the form given by eq. (2/7) in terms of the vectors (VvO^ (which are related to 



the eigenvectors of M by eq. ( |2.12 )) such that au and bk of eq. ( |2. 18|) , (|2. 17j ) satisfy eq.( |2.14 ) or 



equivalently that the T k of eq. ( |2.7| ) are Clifford operators respectively, i.e. satisfy eq.( |2.6D , which 
was assumed before deriving the eigenvalue equation. The corresponding orthogonality relations for 
the eigenvectors of M which are equivalent to the anti-commutation relations for the operators dfc, bk 
lead to further relations between the eigenvectors (see eqs.( [2.4d ) and ( p.42[ )) for special choices of 



the boundary parameters. They simplify the computation of correlation functions and are used for 
projecting to the (++)-sector of B.\ ong . This will be the subject of the sections 9,10 and 11. 

Let us first look at the case where the Hamiltonian is hermitian, i.e. a_ = a* + , (3- = [3+ and 
Q-ziftz £ R» This implies that M has only purely imaginary entries and is also hermitian. So its 
eigenvectors can be chosen to form an orthogonal basis with respect to the standard scalar product. 
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Because M* = — M we have 

which can be directly seen by taking the complex conjugate of the equation M(f>t = A^i". Thus 
after an appropriate normalization of the eigenvectors the orthogonality condition for the eigenbasis 
is equivalent to the relations which are necessary and sufficient to define a set of fermionic operators 



L+l 

3=0 7 

L+l L+l 

E E(tf)7(tf )J = E EWJTOJ = o • (2.2i) 

j=0 7 3=0 7 



Note that due to eq.( 2.19| ) (p k and tfit can always be normalized so that bk and are mutually 



adjoint. For any set of constants Ck G C,c& / the vectors Vfc" = ^( c k4>t + c k l( t ) k) an< ^ 
V'fc = |( c fc^fc ~~ c k 1< ^k) satisfy eqs.( 2.1C| ) and the orthogonality relations 



L+l 

E EWJW)] = C (2-22) 

j=0 7 

and thus the defined in terms of the (V'n)j by eq.(|2.7|) are Clifford operators. If we define \I/ to 
be the (2L + 4) x (2L + 4) matrix consisting of the 2L + 4 vectors i/if, we may rewrite eq.( 2.22| ) as 



vj/tvj/ = 1. This simply reflects the fact that the automorphism group of the Clifford algebra is the 
orthogonal group. Note that here \l/ is not necessary real. However, due to eq.( |2.19| ) \1/ can always 
be made real by tuning the parameters c&. Since ty^ty = 1 implies = 1, we also obtain 

L+l 

E E(€)J(« = Sr (2.23) 

k=0 M 

or in terms of the components of the eigenvectors of M 

EE(«»=2^ • (2.24) 

k=0 M 

Using these equations it is possible to invert eq.( |2.7| ) and eqs . ( 2 . 1 8] ) , ( |2 . 1 7 ) respectively. This is 



necessary to express the spin operators a x ,a y ,a z in terms of ladder operators, which is needed for 
the calculation of correlation functions and of the projection mechanism. We will use this form of 
the orthogonality relations in the next subsection, in order to derive some further relations between 
the eigenvectors. 



12 



If the Hamiltonian is not necessarily hermitian but all of the eigenvalues of M are non-degenerate 
except for the eigenvalue corresponding to the eigenvectors (0, 1, 0, 0, • • •) and (0, 0, • • • , 1, 0), one 
can still show that eqs, (|2.20| ) and ( 2.21 ) remain valid. In general the argument breaks down because 



M is not necessarily diagonalizable. This will become apparent in the sections 3-5. 
Choosing the linear combinations 

0+ = (O,l,O,---,O,i,O) , 0o =(0,l,0,---,0,-i,0) (2.25) 



as the eigenvectors corresponding to the eigenvalue we ensure that they also satisfy the eqs.(2.20) 
and ( gjgjj ). We now check eqs.( |2.2d| ) and fl2.21| ) for the other eigenvectors of M: 



First of all, let (f> + be a right eigenvector corresponding to the eigenvalue A, i.e. 

M(p + = Acp + (2.26) 

Because of M = —M t this eigenvector is also a left eigenvector corresponding to the eigenvalue —A, 
i.e. 

4> +t M = -A(j) +t , (2.27) 

This implies the existence of a right eigenvector cp~ corresponding to the eigenvalue —A. 

Now let (pk and (j>i be eigenvectors corresponding to eigenvalues A^ and A; where we do not 
restrict the real parts of A& and A; to be positive or negative. Using eqs.( 2.26| ) and ( |2.27| ) we get 



4>k l 4>i = ^4>k4>i ■ (2.28) 

A; 



So all products of the form (j)k t (/>i are zero if A& 7^ —A;. This gives eq.( 2.21 ). 

To proof the validity of eq.( |2.20| ) we additionally have to show that in the case — A& = A; the 
product 4>k t< t'i cannot vanish. This can be done by considering <pi*<f>i which is always different from 
zero if <pi 7^ 0. Now due to the assumption of non degenerate eigenvalues the eigenvectors form a 
basis and thus (f>i* can be expressed in terms of eigenvectors (fij. Using eq.( [2.28| ) with — A& = A; we 
have 

O/0f0, = #Vj = E a i#^ = a *40l • (2.29) 
j 

The only term left in the expansion of the product 4>^4>i IS proportional to the product (j>k<j>i due to 
eq.( p.28 ) and therefore cannot vanish. So we can normalize the eigenvectors appropriately in order 



to satisfy equation (|2.20| ). If there are degeneracies in the spectrum of M, the above proof can be 
generalized by using the biorthogonality of left and right eigenvectors. 
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We want to point out that the ladder operators and are not the adjoints of each other in 
general because the relation ( 2.1 9| ) is not valid in general. However, as already mentioned in section 
2.1 it is always possible to perform a similarity transformation in order to achieve b\ = a^. This 
can be seen by choosing an arbitrary real symmetric and orthogonal matrix Vl/'. The transformed 
vectors 



(2.30) 



define a new set of Clifford operators T' k which are now hermitian. Hence the operators b' k = 
\{T' k + + iT' k ~) and a' k = \{T' k + — iT' k ~) form a set of fermionic ladder operators satisfying b k = a' k . 
Since the vectors cf)' k = V4 + i -F V4 _ are no l° n g er eigenvectors of M but of the transformed matrix 

M' = V^M**'* (2.31) 

the transformation ( 2.30| ) corresponds to a similarity transformation of the Hamiltonian Hi ong . 



2.3. Special properties of eigenvectors 

In some cases there are further relations between the eigenvectors in addition to the ones of eqs.(2.20) 



and ( 2.21 ). They are used in the calculation of correlation functions and are even necessary for the 
projection method. First notice that if a z = (3 Z = 0, the matrices M and M 2 respectively take the 
form 

/ n .,. n .,. \ / * n * n . . . \ 



M 



* * 

* * 
* * 

* * 



V 



M 



J 



* * 
* * 

* * 
* o * 



V 



(2.32) 



/ 



where * and both denote 2x2 matrices. Note that * is the notation for an arbitrary 2x2 matrix and 
is not necessarily different from zero. Looking at eq. ( p. 32 ) we see that we can choose an eigenvector 
ijit of M 2 with eigenvalue A|, i.e. 

M 2 ^ = (2.33) 

which satisfies 

$+)f = for i odd . (2.34) 

Now we define ip k by 

= -iA k tP k (2.35) 
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which is also an eigenvector of M 2 with eigenvalue A|. Note that this definition does not work if 
Afc = and thus we have to exclude in the following k = which labels the eigenvectors corresponding 



to the spurious zero mode. Using eqs,( |2.32|) and ( 2.34 ) we obtain 
(^k )i = f° r 1 even • 



(2.36) 



Due to eqs.( |2.33 ) and (2.35) tp k and tp k also satisfy eq.( |2.10| ) and thus we obtain eigenvectors 



v^fc _ V'fc T iV'fc of M satisfying 

(#)? = (-l)*(fr)J • (2.37) 
Therefore for each pair of vectors 4>t, (j)7 satisfying eq.(2.20), there exists a constant c G C such that 



(2.38) 



Now eq.( 2.24 ) can be rewritten as 

L+l 

(i + (-i) i+i )E(^)J(« 

fc=0 



2^ 



and we end up with 



zW)j(« 

fc=0 



for i + J even. 



(2.39) 



(2.40) 



Since eq.(2.38) is not valid for k = and odd L (see eq.(2.25)) we exclude i = 0,u = +, j = L + l, 7 = 
- and j = 0, 7 = + 5 i = L + 1, v = - in the eqs. (p^39|) and (^40|) . 

If diagonal boundary terms are included, and if a_ = a + and /?_ = /?+, the eigenvectors again 
have a special property. In this case M and M 2 have also the form of eq.( |2.32| ), but now * and 
just denote complex numbers and, in place of eq. ( |2,3Sj ), we obtain 



±c'( 



which gives 



L+l 

£( 

fc=0 



(2.41) 



(2.42) 



Both of the equations ( 2.4C| ) and ( 2.42j ) will be used in the sections 10 and 11. 

Note that the proof of eqs.(2.40) and ( 2,42j ) shown above is not valid if there are degeneracies or 
zero modes on top of the spurious zero mode in the spectrum of M. However, one can show that it 
is always possible to build appropriate linear combinations of the eigenvectors corrsponding to the 
same eigenvalue such that eqs.(2.40) and ( |2.42| ) remain valid in addition to eqs. fl2~2Cl ) and ( |2T2T|) . 
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This is not automatically true and therefore in explicit calculations one should take care in choosing 
the right linear combination of eigenvectors corresponding to the degenerate eigenvalues. 

If both conditions a_ = a+,/3_ = /?+ and a z = (3 Z = are satisfied at the same time, then 
both of the equations ( 2.40| ) and ( [2.42 ) can be satisfied simultaneously. By comparing eqs.(2.38) 



and (|2ll| ) we obtain the following relation 



, ^^-{(bu)^ for j even 

(<t>k)t = \ 1 ■ ( 2 - 43 ) 

{ T^(h)f fori odd 

Thus — is either 1 or —1 because otherwise the vectors (jr^ would vanish. 

Let us briefly summarize what we have obtained so far. If H is hermitian, the fact that M = —M l 
leads to eqs.( |2.20| ),( |2,21| ) which are equivalent to the anticommutation relations ( p. 14 ) of the ladder 



operators given by eqs.( |2.18| ), ( 2.17] ). These equations are still valid if H is non- hermitian, but 



diagonalizable. Some additional properties of the eigenvectors have been derived for special choices 
of the boundary parameters (eq. ( 2.4C| ) for a z = j3 z = and eq. ( 2.42| ) for a_ = q + ,/3_ = j3 + ) 



which will be used in the sections 10 and 11. 

3. Calculation of the eigenvectors of the matrix M 

In the previous section we have shown some general properties of the eigenvectors without computing 
them explicitly. This computation is the subject of this section. We will also show how to obtain 
the characteristic polynomial which gives the eigenvalues of M. This polynomial will be treated 
extensively in the following sections. Since M is non hermitian in general we will also discuss the 
diagonalizability of M. 

The eigenvalue problem given by eq.( |2.13 ) is equivalent to a set of recurrence relations. Using 



the notation given by eqs.( 2.1 and (|2,9|) for the eigenvectors (4>^) of M let us first look at the bulk 
part : 



J((^fctf + (^)/ + 2)=±Afe(^)7+l 

-i((^)7 + (^)7 + 2) = ±Afc(^)/ + i 



[l<3<L-2) . (3.1) 



These bulk equations ( |3.1| ) can be decoupled by defining 



<Pi = (<t>i)J + K<i>i)t . ^ = (4)7-K^)t (3-2) 



which gives 
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(3.3) 



Here A = ±Aj. and the functions ipj and Tpj refer to (<pf^)j and {4>^)~j for A = and to (4>k)j an d 
(^fc )j~ ^ or ^ = — -^-fc- From now on we will keep A; fixed and omit all subscripts referring to k. 
Next we treat the left boundary, one obtains 

_ (3.4) 



From the right boundary one gets 



Note that we have excluded explicitly the eigenvectors (0, 1, 0, ...) and (..., 0, 1, 0), which always exist, 
from the set of solutions of the boundary equations above by setting (fo = Tp and Tp L+1 = —Tp L+l . 
Thus we will obtain at most 2L + 2 linearly independent solutions instead of 2L + 4. 



The general solution of the bulk equations (3J3) for A ^ ±i is given by 



ifj = ax 3 + bx 3 , = g(—x) 3 + f(—x) 3 , (3-8) 

where 1 < j < L and up to now a, b,g, f are free parameters which are independent of j. The new 
variable x is related to the eigenvalue A via 

A = ^(x + x- 1 ) . (3.9) 

For A = ±^ the general solution is 

^ = a (±iy + b(±i) 3 j , tp j= g( T iy + f( T iyj . (3.10) 

The four parameters a, b, g, f, the undetermined components (po,Tp , (Pl+i,Pl+i an d the eigenvalues 
A are all fixed by the boundary equations (3.4-3.7) — up to the normalization constants of the 
eigenvectors. Namely, plugging eq. fl3.8| ) into eqs.(3.4)-(3.7) we obtain a homogeneous system of 
eight linear equations with the unknowns a,b,g,f and (po,^, (fL+i^L+i- The condition for the 
existence of nontrivial solutions of this system is given by the vanishing of the determinant of the 
corresponding 8x8 matrix. This defines a polynomial equation in the variable x which yields all 
eigenvalues A. Note that for x = ±1 the 8x8 system of equations has always the non-trivial solution 
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a = —b, g = —f. This corresponds to the zero vector <pj = Tpj = Vj. To compensate this fact we 
divide the polynomial by (1 — x 2 ) 2 . The treatment of the resulting polynomial equation will be the 
subject of the next section (see eq.(4.2)). 

In the following we will show how to obtain the eigenvectors for A 7^ dbi and A 7^ which may 
be viewed as an alternative way to obtain the secular equation. Substituting (3.4) into (3.5) using 
the identity a = <fi\x~ l — bx~ 2 and g = — v^x" 1 — fx~ 2 (see eq. fl3.8|) ) renders b and / as functions 
of if \ and 7p 1 , i.e. 

,2 



1 



1 — x 



a_a_|_ 
x + x _1 



y/2a z - x 1 1 ip± 



x + x 1 



(3.11) 



/ = - 

x 

Thus a and g are given by 
1 



1 



1-x 2 



a_ai-|_ 
x + x~ l 



+ y/2, 



a: 



x + x 



-<Pi 



a_a< 



x + x 



V2, 



a z - x \tpi 



a] 



X + x 



-<PI 



(3.12) 



(3.13) 



a_a_|_ 



a z -x\(p x 



at 



x + x' 



From the right boundary, by substituting (3.7) into (3.6), we see that furthermore 



<PL-1 + 



<Pl-i 



^-V^-4aW + §^ = 



^ + v^-4aW-§^ = 



(3.14) 

(3.15) 
(3.16) 



Using eqs. (|3.11H3.14| ) and eq. (|3.8| ) in the equations fl3.15|) and ( |3.16| ) we get a linear system of 
equations of the form 

/ 

(3.17) 



^21 ^22 



<Pl 



where the Qjj are the following functions of x and the 6 boundary parameters a± , 0±,a z and (3 Z 
(Note that A is a function of x according to eq. (|3.9| )): 

\2" 



n 



11 



x~ L 


' (0-0+ 


1-x- 2 


V 4A 


x L 


'(P-P+ 


1-x 2 


V 4A 



V2(3 Z -x 
V2(5 z -x 



1 \ / a_a+ 
4A 



a - - x 1 1 + 



a_a_ 
4A 



V2a z -x ) + (-1) 



n L (ft-**-)' 
j 4A 2 



16A 2 



(3.18) 
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12 







x 



21 



I -XT 



X 



1 



X^ 
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1 



4A 
4A 



a_a_(_ 
4A ' 

4A ' 

(«+£- ) 2 
16A 2 

16A 2 



x 



X 



a 



7T + (-l) 



4A 



4A 



+ (-1) 



4A 

L [ a - a + 
4A 



+ V2, 



a.. 



+ V2< 



« 2 



-l 



4A 



4A 



+ 



(3.19) 



\/2. 



a 2 — x 



OS-'- 1 '' 



4A 



+ v / 2/3. 



4A 



+ 



_ «2 

V2a z -x) g- + (-l) L 



4A 



4A 



(3.20) 



+ ("1) J 



4A 



+ V2(3 Z 



4A 



+ V2, 



+ 



a_a_)_ 
4A 



+ y/2, 



a? 



X 



The necessary condition to have non-trivial solutions is 
^11^22 — ^12^21 = 



(3.21) 



(3.22) 



This condition is equivalent to the polynomial equation which is obtained from the homogeneous 
8x8 system of linear equations mentioned above. 



The construction of eigenvectors shown here is not valid for A 
show that the eigenvectors for A = ±i can be obtained by 



±i and A = 0. But one can 



l p j = lim (oa^+fr-') , VmM-xY + f{-x)~>) , (3.23) 

x— »±1 a;— >±1 

where a,b,g,f are given by eqs.(3.11)-( |3.14 ). Using de L'Hospital's rule one recovers the form of 
eq. fl3.10l) . The vector components tpi and Tp l are again given as solutions of the 2x2 system ( 3.17 ) 
using x = dbl. 



The solution of the 2x2 linear system (3.17) is straightforward for a given set of boundary 
parameters and a given value of x. It cannot be given in a unique form because some of the fa- 
nlight vanish. We will give the explicit form of the eigenvectors for some special choices of boundary 
parameters in |23j where we are going to calculate the expectation values of crj and crj where j 
denotes the position on the lattice. 

If all f2jj vanish the corresponding eigenvalue is at least twofold degenerate and we obtain 
two linearly independent eigenvectors since (fx and Tp 1 can be chosen independently of each other. 
On the other hand, if a zero of the polynomial is twofold degenerate, it is not clear that all Vlij 
vanish. This comes from the fact that the Hamiltonian is in general non-hermitian and might be 
non-diagonalizable. 
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We would like to point out that the appearance of an eigenvalue A / which is more than 
twofold degenerate would prove that M is non-diagonalizable. This is indeed the case for some 
special choices of the boundaries for a given lattice length L. This can be seen by looking at the 
factorizations of the polynomial obtained in section 5. 



Up to now we have shown how to construct the eigenvectors of the matrix M defined by eq.( 2.11 ). 
For A ^ the components of the eigenvectors are given by 

1, i 



v^)7 = f(^+^) , (0fctf = ~(¥>i-?i) . ( 3 - 24 ) 

where tpj and Tfj have the form of eq. (|3.8|) for A ^ ±| and are given by eq, (|3.23|) for A = ±|. The 
parameters a, b, g and / are given by eqs,( |3.lil - |3.14 ). Finally A has to be determined from eq. ( 3.22 ) 



as a function of x (see eq. (|3.9|)), The vanishing of the determinant in eq.( 3.22| ) leads to a complex 
polynomial of degree 4L + 4 in the variable x which has to be zero. This polynomial will be the 
subject of the next section. 

To obtain eigenvectors corresponding to A = one has to solve the boundary eqs.(3.4-3.7) using 
the bulk solution given by eq,([3.§|) with x = ±i. The calculation is given in the appendix. It turns 
out that besides the eigenvectors (0, 1,0,.. .) and (. . . , 0, 1, 0) of M, which are always present, one 
may have additional eigenvectors for A = 0. Running through the calculation it turns out that this 
happens if = 0. Under this condition two further linearly independent solutions always 

exist. If all relevant boundary parameters vanish, i. e. a_ = a+ = /3_ = /3+ = and if at the same 
time a z = —f3 z for L odd or a z = ^g- for L even we have four additional solutions. Two of them are 
just (1, 0, 0, . . .) and (. . . , 0, 0, 1). Note that the degeneracy of the eigenvalue A = might be higher 
than the number of linearly independent eigenvectors since M might be non-diagonalizable. This 
will be discussed in the appendix by considering the polynomial equation which is given in the next 
section. 

The calculations we have done so far enable us to give a complete set of conditions under which 
M is non-diagonalizable. This is always the case if the degeneracy of an eigenvalue is higher than 
the number of linearly independent eigenvectors. The conditions for the eigenvalue A = are derived 



in the appendix, whereas the conditions for the eigenvalues A 7^ are obtained from eq.( 3.17 ). M is 
non-diagonalizable, if one of the following conditions is satisfied : 

(i) M has an eigenvalue A 7^ which is more than twofold degenerate. 

(ii) M has an eigenvalue A 7^ which is twofold degenerate, but at least one of the is different 
from zero. 
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(iii) A = is an eigenvalue of M, but it is more than sixfold degenerate. 

(iv) A = is a sixfold degenerate eigenvalue of M, but one of the parameters a + , a_ , (3 + , /3_ is 
different from zero. 

(v) A = is a sixfold degenerate eigenvalue of M, but a z ^ —f3 z for L odd or a z ^ ^ or ^ even 
respectively. 

If none of these conditions is satisfied, M is diagonalizable. 
4. The polynomial equation 

Now we turn to the polynomial equation which determines the eigenvalues of M. For later 
convenience, we define a new variable z = x 2 . As can be seen directly from M (eq. ( |2.11| ), one 
of the eigenvalues is always zero, the others are obtained from eq. (|3.9| ) as 

A„ = -g{yfz^ + \Jl/Zn) (4.1) 

Since we have defined the A n to have a positive real part, we will — here and in the following — 
always by definition take the square root which has positive real part. The z n are the solutions of 
the following polynomial which has been obtained from eq. ( p. 22 ) 



p( z ) = _ A {z 2L + 3 + z) + {B + E 2 ) {z 2L+2 + ^ 

(z-iyi 

+ (D + 2E 2 ) (z 2L+1 + z 3 ) + E 2 {z 2L + z 4 ) - 2E {z L+A + z L ) 

+ (—1 + A — B — D) — {-l) L C - 2E 2 ) (z L+3 + z L+l ) (4.2) 

+ C-1 + A-B-D + 2(-l) L C + AE - 4E 2 ) z L+2 + 1 
1 



(z-1) 2 
The coefficients are given by 



q(z) = 



A = 2(— 1 + a-a + + (3+f3- + a 2 , + 01) 

B = (-1 + 2a_o+)(-l + 2/3+/3_) + 4(-l + a_a+)/3 2 + 4(-l + /3+f3-)u 2 z 

C = (a 2 _l3 2 + + a 2 + p 2 _) (4.3) 

D = 2(-l + 2a_a + )/?f + 2(-l + 2(3 + (3^)a 2 z 

E = 2a z p z 

Note that the polynomial p(z) is already completely determined by five complex parameters 
although we started with six parameters in the original Hamiltonian. This can be explained by the 
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existence of a similarity transformation of the form 

/ 



H' = UHU- 1 with U = Y[ h ® • • ■ Ij-i 

3=1 



T j+ i <8> • • • ® Jl (4.4) 



V 



i o 

e 

containing one free parameter e. Here i, stands for the identity matrix at the site j. By using 
this similarity transformation, the four boundary parameters a_,a+,/3_ and /?+ are transformed as 
follows 

a_^ea_; a + — >l/ea + ; /3_ — >• e (5 + — ► 1 /e /3+ (4-5) 

and by choosing a particular value of e, one can always fix one of the boundary parameters. 

The polynomial q{z) has a very special form, because in comparison with a general polynomial 
of degree 2L + 4 many of the coefficients are zero. This changes of course when it is divided by 
(z-lf. 

Observe that the polynomial p(z) has degree 2L + 2 although the diagonal form of Hi ong given 
by eq. ( |2.16| ) has only L + 2 fermionic excitations A n . The reason therefore is the quadratic relation 
between z and A. Since with z also - is a solution of the polynomial, one gets each value of A twice. 
Taking half of them and adding the additional eigenvalue which was already mentioned in section 
2 and explicitly excluded from the set of solutions in section 3 gives exactly the L + 2 fermionic 
excitations. 

Special solutions of this polynomial will be studied in the next section. 
5. Factorization of the polynomial in cyclotomic polynomials 

The study of the factorization properties of the polynomial given by eq. (4.2) represents a very 
interesting mathematical problem. Furthermore, factorizations of the polynomial into cyclotomic 
polynomials which we are going to present below are very important because they allow to calculate 
the whole spectrum and other properties of H\ ong analytically. 

For some special choices of the parameters A, B, C, D and E the polynomial factorizes exactly 
into cyclotomic polynomials. These factorizations were found using the following algorithm: A 

factorized polynomial of degree 2L + 4 of the form 

k 



i=l 



(l-p k+1 z 2L+4 -Z n >) (5.1) 



is expanded for a fixed value of k (corresponding to a fixed number of factors), a fixed value of L and 
for all possible combinations of the rtj. The coefficients of the expanded polynomial f{z) which are 
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functions of the p^, i = 1, . . . , k + 1 are compared to the coefficients of the original polynomial q(z). 
In this way, one obtains a set of 2L + 4 coupled equations for the pi and the coefficients A,B,C,D 
and E which has been solved using Maple. Typically values of L = 5,6,7 were used; for smaller 
values of L the equations do not reflect the general situation because some exponents coincide. For 
larger L however, the number of partitions of 2L + 4 into the rii,i = 1, . . . , k gets too large. Among 
the solutions only those were kept which are valid for arbitrary L and not only for the special L used 
in the calculation, k has been varied between 1 and 4 (i. e. factorizations in 2, 3, . . . , 5 factors were 
studied). For k = 5 (six factors) the program did not run properly - it needed too much memory. 
However, q(z) cannot factorize into a a larger number of factors of the above form Q5.1| ) than six 
with the condition that the corresponding rii appear explicitly as exponents in q(z). In this case the 
only possible combination for the rii would be 

rii = n-2 = L ; ri3 = n.4 = 715 = n§ = 1 (5-2) 

In table I all factorizations which were calculated as described above are listed. The 
factorizations in six factors (entries 14 - 16) were found by solving the system of coupled equations 



for the choice (5^2) of the nj and various choices of the pi. Therefore the list might not be exhaustive 
for the factorizations into six factors. 

The entries 10 to 16 each furnish a one-parameter family of solutions (where the parameter is 
called s) for which a factorization in cyclotomic polynomials appears. Notice however that E is the 
only parameter whose value is always fixed, independently of s. Since E = a z (5 z this means that 
in all cases presented in table I the product of the coefficients in front of the diagonal boundary 
terms is always fixed. Moreover, the entries 10 to 16 provide examples where some of the zeros of 
the polynomial q(z) are always independent of L (e. g. z = s or z = 1/s). We will come back to this 
point in section 7 and in the discussion (section 14). 

Looking at the entries 10 to 16 the remark we made at the end section 3 that in special cases 
the Hamiltonian might not be diagonalizable becomes clear: It is possible to choose the parameter 
s as a function of L in such a way that the polynomial has zeros which are more than twofold 
degenerate. Take for example the case 13 and choose s equal to one solution of 1 — z L = 0. Then the 
corresponding zero of q(z) is threefold degenerate for the value of L chosen above. In this case, one 
can not find more than two independent eigenvectors for M (cf. eq. ( 3.17] )) for the corresponding 
degenerate eigenvalue. Thus in these special cases (which can be constructed analogously for the 
other entries 10 to 16) the matrix M is not diagonalizable. 

For all examples in table I it is possible to calculate the spectrum and the ground state energy 
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Table I. Cases where the polynomial factorizes into cyclotomic polynomials 
A B (-l) i C' D E q(z) 



1.) 




1 


















(1- 2 )(1-2 2L+3 ) 


2-) 







-1 















(1 - z 2 )(l - z 2L + 2 ) 


3.) 




-1 


















(l + z)(l- z L+1 )(l- z L+2 ) 


4.) 




1 


-1 









1 





(l + z)(l-z) 2 (l + z 2L+1 ) 


5.) 












1 

2 










(1 - z L+1 )(l - z L+3 ) 


6.) 












1 

2 










(l-z L + 2 ) 2 


7.) 




1 







1 










(1 - z)(l - z L+1 )(l + z L+2 ) 


8.) 




1 







-1 










(1 - z)(l + z L+1 )(l - z L+2 ) 


9.) 




2 


1 




s + 1/s 










(1 - z) 2 (l - sz L+1 )(l - 1/s z L+1 ) 


10.) 
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exactly. In the following sections we will give the explicit expressions for the ground state energies 
and for some excitations. With the insight gained from these exactly solvable cases we will later also 



treat the general case in the limit of large L [23|. 



6. Examples of exact calculations of spectra for the finite lattice 

Let us now present two examples how to calculate the spectrum of Hi ong from the factorized form 
of the polynomial. We will first take case 4 from table I. The factorized form of the polynomial q{z) 



as 



given in table I suggests the ansatz z = exp(i7r + 2I+1 ) ■ This leads to the solutions 

(j) = nir, n = 1, . . . L (6-1) 

The factor (1 + z) leads to the additional solution z = —1 which means A = 0. The factor (1 — z) 2 
has to be divided out because the fermionic eigenvalues are given by the zeros of p{z) and not of the 
ones of q(z) which is given in table I (also cf. eq. (4.2)). 
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Then the energies of the fermionic excitations are given by A = (twice) and, using eq. ( [4 . 1| ) , 

by 

1 / mr \ 

A n = -sin — n = l,...L (6.2) 



2 \2L + l, 

The energies of the fermionic excitations in the other cases from table I have a similar form. 
However, in the cases 10, 11 and 13 to 16 there is always one solution with roots z = s^ 1 leading 
to a fermionic energy A = ^(v^ + which is — in contrast to the fermionic energies obtained 

in eq. fl6.2| ) — independent of the lattice length L. We will see later [23 that this energy can 
be connected to a boundary bound state. The nature of this state will be elucidated by studying 
the corresponding spatial profiles and by comparing some spatial profiles for special choices of the 
parameters to a Bethe-Ansatz solution of the XX-chain with only diagonal boundary terms [fill. 



This will be described in detail in [23]. 

We now consider the case 9 of table I which is special because the roots of the polynomial all 
depend on a parameter s and in general do not lie on the unit circle. Therefore we will briefly present 
the corresponding results here. For L even we obtain the solutions A = and 

1 /(2n + l)vr ilns \ L 
A n = - sin —, — H — ; n = 0. 1, ... — 

2 V 2(L + 1) 2(L + 1)/ 2 

1 /(271 + lW ilns \ L-2 , . 

An = -sin ■ , - — ■ - n = 0,l,...— - (6.3) 



2 \2(L + 1) 2(L + 1 

for |Jm(ln s)\ < ir. 

For L odd we have accordingly A = 0, and 

1 / nir ilns \ L + 1 

A n = - sin + — — n ' 



2 V^ + l 2(L + 1)/ '"' 2 

1 / nir i In s \ L — 1 , 
A n = -sin — n = 0,l,... (6.4) 



2 V^ + l 2(L + 1 

for < Jm(lns) < ir. For — ir < Im(lns) < one has to interchange the limits of n in the two sets 
of eigenvalues given by eq. (6.4). 

In this example, we can see explicitly how the parameter s appears in the spectrum. The 
argument of the sine is shifted by the s-dependent term • 

For the examples given in table I, it is also possible to calculate the ground state energies exactly. 
The corresponding expressions will be given in the next section. 
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7. Exact expressions for the ground state energies of H t 



long 



Let us first consider the case 4 of table I again. The ground state energy is given by summing up 
all negative eigenvalues of M (cf. eq. (|2.16|) ), this leads to 

E = — > sin = — cot (7.1 

n=l 

In the case 9, also discussed in section 6, the ground state energy for L even is given by 



(L-2)/2 
Z n=0 



(2n + l)7r ilns \ / (2n + lW ilns 

sin I ~t~^ + ^tt tt + sin ' 



2(L + 1) 2(L + 1)/ V 2(L + 1) 2(L + 1 



1 , / Ins \ lcosh 
cosh ' 

2 \2L- 

For L odd we obtain 



icosh = -i wo " 2 ^+ 2 (7.2) 

2 \2L + 2J 2 shi' 77 



2L+2 





2 



7r ln(s) . . ln(s) 
cot cosh h i sinn ■ 



2L + 2 2L + 2 2L + 2 



i r osh ln(s)+i7r 

1 CObll 2L+2_ ^ ^ 



sm ■ 



2L+2 

This expression is indeed real if the original Hamiltonian is chosen to be hermitian. This can 
be seen by imposing hermitian boundary terms in the Hamiltonian a_ = a^_,/3+ = /?* and solving 
the system of equations 

A = 2 = 2(-l + |a_| 2 + |/3+| 2 ), 

B = 1 = (l-2|a_| 2 )(l-2|/5+| 2 ), (7.4) 
(-l) L C = s + \ = 2|a_| 2 |/?+| 2 cos2(£ a + ^) 
where £ Q ,£/3 are the phases of a_ and f3 + , respectively. These equations have only the solution 
|a_| = \j3 + \ = 1, s = exp (±2i(£ Q + ^,g)). For this choice of s the ground state energies given by eqs. 
(7.2) and (7.3) are always real. 

In table II, the expressions for the ground state energies for all cases from table I are listed. 
Again, they correspond to Hi ong . In section 12, we will describe how the corresponding ground state 
energies for the original Hamiltonian H are obtained from the ones for Hi ong . It is remarkable that 
despite the non-diagonal boundary terms all expressions for the ground state energies of Hi ong are 
given in terms of trigonometric functions. The only exception is case 10 where hyperbolic functions 
appear in the expression for the ground state energy (see eqs. (7.2) and (7.3)). However, even in 
this case the model is integrable. 

The virtue of table II is that one can explicitly see how the ground state energy is changing 
with different boundary parameters (we refer to table III of section 12 for some examples in which 
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boundary parameters of the Hamiltonian correspond to a given choice of the A,B,C,D and E). 



This is especially interesting when studying the thermodynamic limit as we will do in [24]. There we 
will show that the Hamiltonian with arbitrary boundary terms corresponds to a conformal invariant 
theory. In particular, if one expands the expressions of table II in powers of one can already 
see that they have the typical form of the ground state expansion corresponding to a conformally 
invariant theory and one can directly read off the conformal charge and the surface free energy for 
the different boundary parameters. 

Notice that the table includes a well-known spin chain: the XX-chain with open boundaries is 
case 10 with s = — 1. 

The s-dependent cases 9-16 are of special interest because the s-dependence is still manifest 
in the expressions of the ground state energy and one can see how this family of ground states varies 
with s. Note that in all cases 10 - 16 the s-dependent terms appear additively to the L-dependent 
part of the ground state energy, (and therewith contributes additively to the surface free energy 
in the expansion for L — ► oo). The physical consequences of the L-independent solutions will be 



discussed in the next paper in connection with boundary states [23]. 

In case 9 however, there is no such simple structure and the parameter s is coupled with L as 
argument of the cosh-term. Indeed the case 9 is special as can be seen already from table I (the 
s-dependence appears in the factors of the polynomial which have degree L + 1, and not in the 
factors of degree 1 or 2 as in all other s-dependent cases). This case will be discussed in detail in 



the following two articles [23 



8. Example of a spectrum of Hi ong with asymmetric bulk terms 

Up to now, we did not look explicitly at Hamiltonians with asymmetric bulk terms which we already 



mentioned in the introduction. Recall that we can map such a Hamiltonian H given by eqs. (1.2) to 
a Hamiltonian H of the form given by eq. ([hi]) which has symmetric bulk terms and L-dependent 
boundary terms of the special form given by eq. (1.3). Using the methods previously described we 



can solve the eigenvalue problem for the Hamiltonian Hi ong and in this way obtain the spectrum 
of Hiong where Hi ong is obtained from H in the same way as Hiong from H by adding one site at 
each end of the chain. We will carry this out for one example. We choose a Hamiltonian Hi on g 
whose transformed Hi ong has boundary terms corresponding to the case 9. Notice that this is the 
only factorizable case where this can be done independently of L if p ^ q, i. e. the analysis can 
simultanously be carried out for all Hamiltonians H of the chosen type with arbitrary length. In all 
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Table II. Ground state energy of Hi ong for the cases from table I 
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1 rot - 

2 LUU 2L+2 

2iI^-|( S 1/4 + ^ 1/4 )-i( S 1/4 - S - 1/4 ) 
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other cases starting from one factorizable case and changing the length L would result in boundary 
parameters belonging to another (perhaps not even factorizable) case. 

One possible choice for the boundary parameters in the case 9 is given by (as can be directly 
seen from eq. (4.3)): 



a-a+ = l; 13+13- = 1; {-l) L (a 2 _ /3 2 + + a\0t) = s + a z = = (3 Z 



1 



Expressing a_ in terms of a+ and /?_ in terms of /3+ and using eq. ( |1.3| ), we obtain s = 
(-1) L (j^J Q 2 - 2L (or 1/s = (-l) L Q 2 ~ 2L ). Using the results obtained in eq. (6.3) we get 



for the fermionic excitations for even L 



. 1 . I (2n + 1)it . 

A = - sin — ± 1 

2 \ 2L + 2 



21nQ + ln^ 



L + l 



InQ 



n = 0, . . . 



L — 2 



2* 



A = — cosh 
2 



21nQ + ln 
L + l 



lnQ 



The quasi-momenta ^l+I^ are shifted by the constant i 



21nQ+ln 



InQ 



U) 



The second L- 



independent term of this constant is typical for asymmetric bulk terms as will be seen in [24]. A 
similar expression has been obtained in || for the fermionic excitations of H with totally asymmetric 
bulk terms (p = 1, q = 0), a_ / 0, /3+ / and a + = (3- = a z = [3 Z = 0. 

The ground state energy is given by summing up all negative eigenvalues of M, leading to 



in 



Eq 



cosh(- 



L+l 



InQ) 



2sm(^) 



i-2) 



Observe that for | = 1 we obtain the previous expression calculated for the case 9. 
9. Projection method and the a x 1-point functions 

Up to now we have dealt only with the Hamiltonian Hi ong given by eq. (2.1) which was obtained 
from H (see eq. (1.2)) by adding one lattice site at each end of the chain. In this section we will 
explain how the spectrum of H is related to the spectrum of Hi ong . 

Since, as mentioned in section 2, ctq and cr£ +1 commute with U[ ong the spectrum decomposes 

into four sectors (++, H — , — h, ) corresponding to the eigenvalues ±1 of erg and cr£ +1 . The 

eigenvalues and eigenvectors of H are related to the (++)-sector in the following way: 

If \E) is an eigenvector of H corresponding to an eigenvalue E then \Ei ong ) = |+) ® \E) ® |+) is 
an eigenvector of Hi ong corresponding to the same eigenvalue E, where 0"ol+) = |+) and \Ei ong ) is 
element of the space C 2 ®C <S>C 2 . Thus the whole spectrum of H is contained in the spectrum of 
Hiong projected onto the (++)-sector. Since the dimension of the of (++)-sector is 2 L we conclude 
that the spectrum of H is identical to the spectrum of Hi ong projected onto the (++)-sector. 

Before describing how we will proceed to project to the (++)-sector, we make some definitions 
which are needed later. First, we want to remind the reader that M defined by eq.( p.ll ) always 
has a twofold degenerate eigenvalue corresponding to eigenvectors (0, 1, 0, 0, • • •) and (0, 0, • • • , 1, 0) 



of M. Using the Clifford operators rj given by eq. (|2.3|) , we now define the corresponding ladder 
operators 



(4 



IT 



L+l 



)/2 



«o = (tq +ir L+1 )/2 



(9.1) 
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We also define the vacuum representation with the lowest weight vector \vac) by 

a k \vac) =0 Vfc . (9.2) 

Because we are interested in eigenstates of (Jq and cr£ +1 we define the vectors 

\v ± ) = ±={\vac)±\Q)) , (9.3) 

where |0) = bo\vac). Observe that 

a^} = ±\v ± ) . (9.4) 

Now we will proceed in three steps. In the first step (subsection 9.1) we will show using some 
algebraic considerations that the vectors \v^} are also eigenvectors of cr£ +1 . It will turn out that the 
eigenvalues of <r£ +1 corresponding to the eigenvectors \v + ) and \v~) always have opposite signs, i.e. 

al^) = ±r l \v ± ) (9.5) 

with rf = 1. The value of rj plays a crucial role in the following. We will also show that the 
(++)-sector consists either of the states 

r 

Y\_ bkj\v + ) with r even (9-6) 
J'=l 

or 

r 

Y[b k .\v~) with r odd , (9.7) 

i=i 

where < kj < Since kj ^ 0, the creation operator of the spurious zero mode defined by 



eq.(|9.1|) does not appear in (9J3) and (9^7). The groundstate of H corresponds to \v + ) or to bk\v~), 
where b^ denotes the creation operator corresponding to the fermionic energy with smallest real part 
and k ^ 0. 

Whether the (++)-sector consists of the states (|9.6[ ) or (^) will be shown to depend on the 
value of r\ (see eq.( |9.5[ )) which will be calculated by computing the expectation value (v + \a'f j+1 \v + ). 
Note that we define (v \ via 

(v^l = — j=({vac\ db (vac\ao) , (9-8) 
v2 

where (vac\ denotes the left vacuum of Hi on g, i.e. 

(vac\b k = . (9.9) 

Note that if Hi ong is not hermitian, (vac\ is not equal to |uac)' in general. 
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In the second step (section 10), we will show how to calculate the one-point function 

f(J) = (v + Wj\v + ) (9.10) 

In the present context, this is done merely for technical reasons in order to calculate /(L+l), however 
it will be essential in another context. Namely, the explicit calculation of f(J) will be presented in 



the following paper [23|. We would already like to remark that the calculation of f(j) is similar to 



the calculation of the two-point function 

g(i,j) = (v + \ f<#\v + ) . (9.11) 



Note that due to eq,( |9.4| ) one already sees that f(j) = g(0,j). We show in section 10 that the 
functions /(j) and g(i,j) are both given by Pfafnans of submatrices of the same matrix. We will 
furthermore generalize these considerations to states of the form ( |9.6| ) and ( |9.7| ). These results also 
apply to H although we started with the larger space of states of H[ ong . 

Using eq.( 2.40| ) and eq. (|2.42| ), we will obtain determinant representations for f(j) and g(i,j) 



which can be treated analytically in the calculation of f{L + 1). This calculation can be found in 
section 11 and will be the third step. 

9.1. Algebraic considerations 

In this subsection, we will show that the (++)-sector consists either of the states given by fl9.6|) or 



d9.7|) and clarify the role of rj given in eq. (|9.5|) . Because 



(^ = &i = ^M0 (9.12) 
which can be seen directly from the matrix M, we obtain the commutation relations 

[°L+i,h} = K+i,ak) = {o-lb k } = {a%,a k } = ,Wk^0 (9.13) 
from eqs. ( gig) and ( [ml) and therefore 



H,N k ] = [o- x L+1 ,N k ] = ,Vfc#0 . (9.14) 
Due to eq.( |9.14| ) the vectors Iv^) have to be eigenvectors of <r£ +1 , i.e. 

ol fl |u ± >=»r b |u ± > • (9-15) 

Thus the sector containing \v + ) respectively \v~) is well defined and given by the value of r] + and 
rj~ respectively. Note that eq.( 9.15| ) is not as precise as eq. (|9.5| ) because eq.( p.5| ) implies rj + = —r]~. 
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Due to eg. (gig) we can make the following statement concerning the vectors of a given sector: 
If an arbitrary vector \v) is element of the (± e)-sector where e S {+, — } then b^v) with k ^ is 
element of the (=F e)-sector. Now one has to distinguish two cases: 

First, if \v + ) is an element of the (++)-sector, i.e. rj + = +1, then all the states given by (|9.6[) 

are also elements of the (++)-sector. The vector \v~) then has to be an element of the ( )-sector 

because otherwise the ( )-sector would be missing in the space of states which is not the case. 

Thus we have rj~ = — 1. 

Second, if \v + ) is an element of the (H — )-sector, i.e. r] + = — 1, then \v~) has to be an element 

of the (— +)-sector. Otherwise \v~) would be an element of the ( )-sector and there would be no 

(++)-sector. As a consequence we have that r]~ = +1 and all the states given by fl9?7| ) are elements 
of the (++)-sector. 

In both cases, the subspace spanned by the vectors ( |9.6[) respectively ( |9.7D has dimension 2 L and 
thus they form a basis of the (++)-sector. Because the values of tj + and rj~ always have opposite 
signs, we will use the variable r\ defined by eq.( |9.5D in the following. 

10. 1- and 2-point functions of a x 

In this section show how to compute the 1- and 2-point functions for the erj operator with respect 
to the states \v^) defined by eq. (|9.3| ) and (v^l defined by eq. (|9.S|) following the way of Lieb, Schultz 
and Mattis (LSM) Q. The computation of the correlators of the operators erj can be done similarly. 
In this paper, however, we are only interested in the value of {v + \a C j\v + ) at the particular point 
j = L + 1 which is calculated in section 11 using the results of this section. Due to eq.( |9.5D we 
thereby obtain the value of rj which is needed for the projection from Hi ong to H as described in 
the previous section. The calculation for general values of j will be part of the second paper of this 



series | 23fl . 



The main difference in comparison to the problem of LSM |l]] is that in our case the 1-point 
functions do not vanish because of the non-diagonal boundary terms we are considering. They can 
be calculated in the same way as the 2-point correlators. At the end of this section, we will briefly 
remark how to compute the 1- and 2-point functions for excited states of Hi ong of the form fl9.6| ) 
and fl9.?1 ) . This is not necessary for the calculation of rj defined in eq. ( |9.5| ) , but it is needed for the 
calculation of the 1- and 2-point functions for the eigenstates of H. Note that if r\ = — 1 the ground 
state of H corresponds to an excited state of Hi ong . 

We now proceed to the calculation of the 1- and 2-point functions. Writing <rj in terms of the 
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T i defined in eq.([T3) we obtain 



| a] | «±) = ±(-i)V | r " J] r+r fc -r+ | V ±) , 

k<j 

which is up to the sign exactly the 2-point function (v^ \ (JqCTj \ v^}. In general 



(10.1) 



( U ± i ^ of i = (-iy-*<«= 



(10.2) 



fc<j 

Of course i denotes y/—l. Using eq.( [2.24| ) the Tj can be expressed in terms of ladder operators 
and bk, i.e. 

L+l 

r i=E(€)^ + (^)^ • ( 10 - 3 ) 

fc=0 

Because (4*0 )j is only different from zero if either j = and /i = + orj = L + l and fJ- = — (compare 
ls.(fn8|),(|2~T7l)) we have 

~~~ (10.4) 



vac) 



eq.Q an d eqs.(^18|),(|2~T7|)) we have 

I rr 11 r+T k T+ | = (vac | r~ [] r+r^r/ 

For simplicity, we will denote (vac\0\vac) by (O) in the following where O denotes an arbitrary 
operator. Using Wick's theorem we are left with the calculation of the Pfaffian of the antisymmetric 
Matrix A, i.e. 



(nil 

k<j 



T k T k y 



PfA 



(10.5) 



where 



1 (t-t? +1 ) (rr T r +1 ) 

( r ii r r) ° ( T t+i T i+i) 

( T i+2 T i ) ( T i+2 r i+l) ( T i+2 T i+l) 





~T? +2 ) • 


•• (n4) \ 




l T i+2) ' 


















(10.6) 



V krfrr) (r/r+i) (r+rr +1 ) ( r +r+ 2 > ••• y 
We want to remind the reader that the Pfaffian of a In x 2n antisymmetric matrix A with matrix 
elements a^- is defined by 

P ^ = ion X! S g n ( <T ) a <r(l)(T(2)a CT (3) ( T(4) • • • a CT (2n-l) CT (2n) j (10-7) 

where S2n denotes the symmetric group of degree In. 

The expectation values of the basic contractions of pairs which form the entries of A are evaluated 

using eq.( 10.3j ). Due to the property of (vac\ (cf. eq.( |9.9[) ) we obtain 

L+l 

W> = £(«W)i ■ ( 10 - 8 ) 



k=0 
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In general, no further simplification is possible. 

Nevertheless there exist two special cases where the calculation of the Pfaffian can be reduced to 
the calculation of a determinant. Namely, this is possible if no diagonal boundary terms are present 
or if a_ = a+ and (3- = /?+. This reduction uses the additional relations for the eigenvectors given 
by eqs.( 2.4C| ) or ( p.42| ) respectively. Note that in this paper the general relation (PfA) 2 = det A is 
of no use because we are exactly interested in the sign of rj. 

In the absence of diagonal boundaries we can use eq.( 2.40| ) to simplify eq.( 10.8| ). Using similar 
arguments as LSM we also obtain a determinant representation for the correlation functions, however, 
the contributing contractions are different from theirs. In fact, using eq.( |2,40 ) results in 

(rf rj) = for i + j even (10.9) 

and thus the correlation functions are given by subdeterminants of the (L + 1) x (L + 1) matrix 

( (^4) (Tor{) (To 4) 

(44) <r+rf> (44) ■■■ 

(^4) ( T 2 T l) ( T 2~4) 



D 



V 



(10.10) 



Denoting by D*- the matrix D after elimination of the first i rows and columns and the last L + l—j 
rows and columns we can write 



PfA = det Dj 



(10.11) 



where 



i ij 



— i if i even and j odd 
i if i odd and j even 
1 otherwise 



(10.12) 



The calculation of PfA also simplifies if a+ = a_ and /3+ = In this case we can utilize eq.(2.42) 
to obtain 



(44) = o . 

This again results in a determinant representation of PfA, i.e. 

(44 + d (n4+2) ■■■ K4) 

■ (4i4) 



PfA 



( T i+\4+l) ( T i+l T i+2) 



( T j-l4+l) ( T j-l4+2) ( T i-1 T /) 



(10.13) 



(10.14) 
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Thus in eq.( |10,5 ) we are left with the calculation of the subdeterminants of a (L + 1) x (L+l) matrix 
G with elements 

GH = (r fc l 1 r+) , (10.15) 

where 1 < k < L + 1 and 1 < I < L + 1. 

By simple modifications one can generalize the results of this section to excited states of the 
form (|9,6| ) or (^^). The argument runs as follows: Any state given by elementary excitations can 
be regarded as the vacuum state \vad) of a new set of ladder operators, where the and of the 
excited fermions are interchanged. This corresponds to an interchange of the eigenvectors 4>7 and 



4>^ (see eqs. fl2.17 ) and ( 2.1§| )). Thus the calculation of correlation functions for the states in fl9.6|) 
and ( |9.7j ) is equivalent to the calculation of the correlation functions for the states \v + ') defined in 
the same way as \v + ) in eq.(fO) but now using 



\vac 



UK 

3=1 



vac) 



(10.16) 



as vacuum state where r is even or odd respectively and k ^ 0. The left vacuum defined by eq. fl9.9| ) 

(10.17) 



has to be modified analogously. As consequence, we only have to replace eq.( 10.8 ) by 



W>= E 

h uncx< 



+ E 



11. Calculation of (<r£ +1 ) 



If the diagonal boundary terms are absent or if a+ = a_ and j3 + = j3- , we can make use of eqs.( |10.lT| ) 



or (10.14) respectively in order to calculate 

v = (v + W x L+1 \v + ) ■ (n-i) 

Recall that we need the value of r] to decide whether the (++)-sector of the space of states of Hi ong 
is given by the states of the form ( |9.6j ) or ( |9.7D . 

Thus if no diagonal boundary terms are present, we have to calculate detD where D is defined 
by eq. ( |10.10| ) . This determinant can be written as the product of two determinants 



(-1) L det D = det A g det A„ 



(11.2) 



where 



(^r)o 
{<t>i)t 



J 2 10 



<>3 )o 
(03)2 



(0 2 + )r 



4)s 



V 



(11.3) 
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The factor (— 1) L is due to permutations of rows and columns which can be seen by comparing 
the product A g A u with D. Without loss of generality we may assume that the eigenvectors are 



A A* 



normalized in such a way that they satisfy eq.( 2.38| ) with c = 1. Using eq.( 2.40 ) we then find 

1 , A U A*=-1 . (11.4) 

Due to the special form of the matrix M (see eq.( |2,32| )) and to eq,( 2.38| ) the two matrices A g and 

via 



A u are related by a matrix M 9 ^ u 

A* 



A- 1 M 5 ^ n A 5 



(11.5) 



where A 



kk> 



A-kSkk' with k > . The matrix M 9 ^ u is given by 
/ 



M 



for I odd or 
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g^u 
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—Ir* 
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(11.6) 



(11.7) 



/ 



(11. 



for I even, where we denoted by G 1 the matrix G of eq.(2.5) with the second row eliminated and 
by K' the matrix K of eq.(2.5) with the first column eliminated. Note that the eliminated rows and 
columns contain only entries which are equal to zero. 

Using the eqs.( ll.Ej ) and ( |11.4j ) in eq.( 11.2j ), we obtain 

detD = (-l) L det A _1 detM g ^ u . 

The value of det M 9 _, M can be computed in an elementary way 

-(a_/9+ + a +/ 3_)/4 L+1 for L odd 
-i(a_ (3 + + a + /5_ ) /4 L+1 for L even 

Plugging this into eq.(ffOT^), we end up with 

L+1 a-f3+ + a + (3- 



detM 



(11.9) 



PfA = (-i) 1 



4^+1 n M o^ 



(11.10) 
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If a_ = a + and (3_ = /?+, we obtain the same result by using equation ( |10,14 ) and performing a 



similar calculation. Combining eqs. (|101) , (|10^) and (|Ta5D with eq.( pTTo| ) we are left in both cases 
with 

,^>! +1 l^(-D"'5§i^ • (ii.il) 

Notice that in these cases the expression of r\ does not depend on a z and (3 Z . It is not possible 



to calculate the product of all eigenvalues in eq. flll.ll ) in general, but the squared product of 



eigenvalues can be calculated from det M' where M 1 denotes the matrix M with the second and the 
last but one row and column eliminated. In both cases this yields 

det A 2 = (-l) L+1 detM' = (a_/3+ + a+/3_) 2 /4 2L+2 . (11.12) 

If the Hamiltonian is hermitian, the product of eigenvalues is simply given by n^o^-fc = \ a -P+ + 
aq-/?_|/4 +1 . This also holds for Hamiltonians with only real entries because in these cases the 
product is real and positive. Thus, we can apply equation ( |1 1 . 11 ) directly. Otherwise one would 
have to know all the eigenvalues explicitly. 

Note that the value of r] can only change by variation of the boundaries if one crosses a point in 
the parameter space at which an additional mode with Re A/% = and k ^ exists. This is due to the 
fact that the eigenvalues of H and H[ ong are continuous functions of the boundary parameters. But 
at a point satisfying the above condition the value of r] is not well defined because the corresponding 
ladder operators at and bk are not well defined as already mentioned in section 2.1. Therefore one 
can aquire a change of sign in r/ by passing through such a point. 

If H is hermitian the condition to have a mode with ReA^ = and k ^ is equivalent to the 
existence of an additional zero mode. The presence of such a zero mode corresponds to the root 
z = —1 in the polynomial given by eq.(4.2). This implies the condition + = 0. Thus if 

H is hermitian, i.e. a+ = a*L,f3+ = (3*_ , we have to only distinguish the two regions Re(a_/3+) > 
and Re(a_/3+) < 0. Thus, we conclude that if H is hermitian we obtain the following expression for 

V 

rj = (-l) L+1 sign(Re(a_^+)) . (11.13) 



The results of this section, namely eqs. Qll.llj ) and ( 11.13 ), allow us to calculate the ground state 



energy of H for the exactly solvable cases of table I. This will be the subject of section 12. 
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12. Ground state energies for the Hamiltonian H in the exactly solvable cases 

In section 9 we have shown that the ground state of H corresponds either to \v + ) or to b[ owes t\v~) 
where bi owes t is the creation operator corresponding to the fermion energy with the smallest real 
part which we will denote by 2Ai owes t in the following, and the \v^) are defined in eq. ( |9,3| ). Which 
of these two states corresponds to the ground state depends on the eigenvalue rj of \v + ) with respect 
to the operation of The eigenvalue rj is either +1 or —1. If i] = 1, then the ground state 

corresponds to \v + ), and the ground state energy of H is equal to the ground state energy of Hi ong . 
If rj = —1, then the ground state corresponds to bi owes t \ v~) and the ground state energy of H is 
given by the sum of the ground state energy of Hi ong and 2Ai owest . For the exactly solvable cases, 
the ground state energies for Hi ong are already contained in table II. 
If at least one of the following conditions is satisfied: 

a) Hi on g is hermitian 

b) Hiong has no a z boundary terms (a z = = (5 Z ) 

c) a_ = a + and /?+ = /?_, 

the value of r\ can be easily calculated by using the explicit formulas of section 11. In the cases b) 



and c), the expression for r\ is given in eq. (11.11) in terms of a_, /3+, /?_ and a+ and the eigenvalues 
Afc of Hi ong . In the case a) where Hi ong is hermitian, r? is given by eq. ( |11.13j ) in terms of a_ and 
/?+ alone. In the other cases, one would have to calculate the Pfaffian of the matrix A given by eq. 
( |10.6| ) using different methods than the ones we used in the sections 10 and 11 to decide which of 
the states \v + ) or bi\v~) corresponds to the ground state of H. 

To determine the ground state energy for a given Hamiltonian H of the type a) — c) one has to 



calculate the expression for 77 given by eq. ( |l 1.11 ) or eq. ( 11.13 ). Analytically, the ground state energy 
for H can be calculated with our methods only for the cases given in table I where the polynomial 
factorizes into cyclotomic polynomials and where one knows the whole spectrum of Hi ong . These 
cases are given in terms of the parameters A, B, C, D and E. In order to put our machinery to 
work, we need the corresponding parameters a_, a+, /?_, /3+, a z and (3 Z . Since the transformation 
from A,B,C,D and E to a_, a+, /3_, /3+, a z and /3 X (which is given in eqs. (4.3)) is non-linear and 
leads from 5 to 6 variables, the choice of the parameters a_ , a + , /3_ , [5 + , a z and f3 z for a given set 
A,B,C,D and E is not unique, creating some freedom of choice. 

For all the cases listed in table I, we solved the eqs. (4.3) for a_, a + , /?_, a z and Z and 
allowed only solutions which additionally satisfy one of the conditions a)-c) above. These solutions 
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are listed in table III. The choices of boundary parameters obtained from the ones given in table III 
by application of an obvious similarity transformation to the Hamiltonian H such as reflecting the 



Hamiltonian in the middle of the chain or applying a transformation of the form given by eq. (4.4) 
to the Hamiltonian are not explicitly listed. 

Table III should be understood as follows: Let us first comment on the choice of the signs in the 
cases where we give two alternative signs for the boundary parameters. The signs in the third and 
in the fourth column can always be chosen independently of each other. However, in the upper half 
of the table, the signs of a_ and a+ cannot be chosen independently (see e. g. case 2 or 4) whereas 
in the lower part of the table, they can be chosen independently (e. g. in the cases 9 even and 9 odd 
as indicated by the condition a + j3- = ±1). 

Now we turn to the one-parameter families depending on the parameter s. Here, the conditions 
a)-c) often lead to restrictions for the value of s which are indicated in the fifth column. The 
choice of a_,a+,/3_ and /?+ given in the fourth column determines a_/?+ + a+/J_. In most of the 
cases from table I the product of all values of is positive (or zero), and consequently the sign 
of a_/? + + ot + (3- multiplied by (— 1) L yields the sign of rj. If rj < 0, the ground state energy of H 
is given by Eo(Hi ong ) + 2Ai owest where Eo(H[ ong ) can be taken from table II. If r) > 0, the ground 
state energy of H is given by Eq(H) = EQ(Hi ong ). However, there are some cases where eigenvalues 
with vanishing real part but non-vanishing imaginary part appear in the spectrum of Hi ong . This 
happens for example in case 10 for negative values of s and may in general happen in the cases 9, 10, 
14 and 15. In these cases, the sign of the product of all eigenvalues is not uniquely defined. Here, 
it is impossible to decide which of the two vectors \v + ) and biowest\ v ~) corresponds to the ground 
state of H. 

In some of the cases, we always find ct-/3+ + a+/3_ = 0. Here, an additional zero mode appears 
in the spectrum of H[ ong as already mentioned at the end of section 3. Therefore the energies of 
\v + ) and b n \v~) where b n is the creation operator for the additional fermion with energy zero are 
the same and the ground state energy of H is again given by Eo(Hi ong ). This is also indicated in 
table III in the last column. 



13. Guide 

In this article, we explained how to diagonalize the XX-quantum spin chain of length L with diagonal 



and non-diagonal boundary terms defined in eq. (1.1). Here we give a resumee of our method which 



the reader may use as a guide on how to use our results. This guide should be seen as a user-friendly 
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Table III. Exactly solvable cases from table I. Details needed for the projection method 
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cooking recipe. It has two parts, the first deals with the spectrum, the second with the eigenvectors. 
As one will notice, the guide does not follow the sections in a chronological way. 

In order to find the eigenvalues and -vectors of H we start by considering a different Hamiltonian 
Hiong which is obtained from H by appending two additional sites and L + 1 (see eq. (1.4)) so 
that the expression for Hi ong is bilinear in Majorana (Clifford) operators, see eq. (|2.4[) . Hi ong can 
be diagonalized in terms of free fermions, fixing the representation we are working in. The spectrum 
and the eigenvectors of H in the Fock representation can be retrieved from the ones found for Hi ong 
by a projection method described below. 

13.1. Eigenvalues of Hi ong 

The diagonalization of Hi ong is described in section 2. The spectrum is given in terms of L + 2 single 
fermionic energies 2A n (see eq. ( [2.16D ). The values of A n can be obtained from a (2L + 4) x (2L + 4) 
matrix M (eq. ( |2.11[ )). Since M = — M', the 2L + 4 eigenvalues of this matrix appear in pairs 
±A n . The necessary L + 2 eigenvalues are taken by convention as the values with positive real part. 
As explained in the text, zero is always an eigenvalue of M. This corresponds to a fermionic zero 
mode. As long as we consider Hi ong , the zero mode Ao = appears in the spectrum and the ground 
state is at least twofold degenerate. As we are going to see, the zero mode does not appear in the 
spectrum of H, therefore we are going to call it spurious zero mode. However, the eigenvectors of 
M corresponding to the spurious zero mode will be needed in the derivation of the eigenvectors of 
H. 

The values of A n can be expressed using eq. (|4.l| ) in terms of the zeros of a complex polynomial 
of degree 2L + 4 (see eq. (4.2)). To find the zeros of the polynomial analytically, we have looked in a 
systematic way for factorizations of the polynomial into cyclotomic polynomials. We determined all 
possible factorizations up to five factors and found some examples for factorizations in six factors. 
These results are listed in table I (the parameters A, B, C, D and E appearing in table I are defined 
by eq. (4.3) in terms of the boundary parameters of the Hamiltonian). For the cases where we did not 
find any factorizations of the polynomial, the zeros of the polynomial and therewith the fermionic 
energies can still be calculated numerically. Since the polynomial has degree 2L + 4, this is much 
easier than a straightforward numerical diagonalization of the Hamiltonian which has dimension 
2 L x 2 L . 

By studying the solutions of the polynomial equation (4.2), we find special L-independent 
solutions in some cases. They correspond to boundary bound states as will be shown in [|^] . 
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The ground state energy of Hi ong (which is by convention the energy with the smallest real 
part) is obtained in eq. ( [2.16D by subtracting the Fermi sea. In table II we listed the corresponding 
expressions for the ground state energies of Hi ong (which are at least twofold degenerate) for the 
cases where the polynomial factorizes into cyclotomic polynomials. Some properties of the ground 
state energies will be discussed in section 14. 

In section 6, we give the expressions for the spectrum of M in some of the "exactly solvable" 
cases. A list of the ground state energies of Hi ong for all "exactly solvable" cases can be found in 
section 7. Section 8 contains the spectrum of M and the ground state energy of one example of a 
Hamiltonian with asymmetric bulk terms which can be treated with the results developed in this 



article by using the similarity transformation between the Hamiltonian given by eq. (1.1) and the 



one given by eq. (1.2). This transformation changes the boundary parameters according to eq. ( |1.3| ). 



13.2. Eigenvalues and ground state energy of H 

Finding the eigenvalues of the original Hamiltonian H is more involved. As shown in section 9, to 
find the spectrum H we have to look at an even or an odd number of fermionic excitations with 
respect to the lowest energy of Hi ong . We disregard the spurious zero mode in the calculation of 
the number of the fermionic excitations. Whether one has an even or an odd number of fermionic 
excitations in the spectrum of H depends on the value of a parameter r\ defined by eq. (|9.5| ) which 
is either +1 or —1 (see section 9 for details). The way it is computed is going to be explained below. 
If 7] = +1, the spectrum of H consists of an even number of fermionic excitations with respect to 
the ground state energy of H[ ong and the ground state energy of H is the same as the one of H\ ong . 
If rj = —1, the eigenvalues of H are given by an odd number of fermionic excitations and the ground 
state energy of H is the sum of the ground state energy of Hiong and the fermionic energy with the 
smallest real part which we call 2Ai oweS f. 

If on top of the spurious zero mode another fermionic excitation is zero, the ground state energy 
of H is non-degenerate and identical to the ground state energy of Hiong and the spectrum of H 
is given by all even and odd combinations of fermionic excitations. If a second fermionic excitation 
is zero, the whole spectrum of M is twofold degenerate. So in these cases one does not need to 
calculate the value of rj. 

At this point we restrict our discussion to the cases, where we have derived explicit formulas for 
the parameter rj: 

a) H^ng (and therewith H) is hermitian 
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b) H[ ong has no o z boundary terms (a z = = f3 z ) 

c) a_ = a + and /?+ =13-. 

For the other cases this guide is not sufficient since they are much more complicated and we have 
not obtained simple formulas for the parameter rj. 

In case a), 77 is given by eq. (|11.13[) ). Notice that only two boundary parameters appear in the 
expression for r\. In the cases b) and c), the expression for rj is given by eq. ( ll.ll| ) in terms of the 
parameters of the non-diagonal boundary terms and the eigenvalues of M. 

13.2.1. Analytical results for the ground state energy of H 

If the Hamiltonian additionally belongs to one of the " factorizable" cases, the ground state energy 
of H can be calculated analytically. These cases are listed in table III. To calculate the ground state 
energy of H for a particular choice of boundary parameters given in the third and fourth column 
of table III, one has to proceed as follows: First one checks the value of 2Ai owes t given in the last 
column. If 2Ai owest = 0, the ground state energy of H is identical to that of Hi ong (which is listed in 
table II). In the cases where 2Ai owes t 7^ 0, one has to know the value of rj to obtain the ground state 
energy of H. This value is obtained by using formula ( 11.11] ) in the cases b) and c) (the eigenvalues 



A„, are given by the zeros of the factorized polynomials listed in table I via eq. (^lj)) and formula 
(11.13) in the case a). 



The fermionic energy with the smallest real part 2Ai owes t which has to be added to the ground 
state energy of Hi ong if rj = —1 is listed in the last column of table III. If 77 = 1 the ground state 
energy of H can be taken directly from table II. 

Many of the exactly solvable cases depend on an arbitrary free parameter s (see table I and III). 
In table III, these s-dependent cases can be separated into two categories. For the cases 11-13 
and for some choices of the parameters in the cases 9, 10, 14 and 15, the conditions a)-c) fix the 
parameter s to some particular value which can be found in column 5 of table III. For the cases 9, 10, 
14, 15 and 16 there are also possible choices of the boundary parameters where this is not the case. 
In the examples 9, 10, 14 and 15 it may happen that one cannot make a definite statement about 
the value of rj, if s is chosen in such a way that the s-dependent eigenvalue of M has a vanishing 
real part, but a non-vanishing imaginary part. The reason lies in the fact that our convention to 
choose the fermionic energies as the ones with positive real part becomes ambiguous in this case. 



13.2.2. Numerical calculation of the ground state energy of H 

Even if the Hamiltonian one is interested in does not belong to one of the factorizable cases, but 
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fulfils conditions a), b) or c), one can still use the formulas ( 11.13] ) and ( 11.11] ) to decide what the 



ground state of H is. If H is hermitian, r\ can be read off directly from eq. ( |11.13| ), in the cases 
b) and c) one additionally needs the spectrum of M to compute the value of r/ (see eq. ([11.11]) ). 
The eigenvalues of M can be calculated numerically by solving the polynomial equation (4.2) or by 



diagonalizing M numerically. Inserting them into eq. ( 2.16| ) yields the ground state energy of Hi ong . 



13.3. Eigenvectors of Hi ong , H and M 

Up to now we have described how to find the eigenvalues and the ground state energies for Hi ong 
and for H. Let us now turn to the eigenvectors. 

The eigenvectors of Hi ong are given in a fermionic Fock representation (compare eq. p. 16 ). 



The eigenvectors of H are given in the same Fock representation, however they all lie either in 
the even or the odd part of the Fock space where we again to not count the spurious zero mode. If 
the value of n is +1, the ground state of H corresponds to \v + ) (see eq. fl^3])) and all excited states 
are of the form given in eq. (|9.6]). If rj = —1, the ground state of H corresponds to bi owes t\v~) where 
\v~) is also defined in eq. fl9.3|) and bi owest is the creation operator corresponding to the fermion 
energy with the smallest real part. For the exactly solvable cases it can be read off table III. The 



excited states are described by eq. (9.7). 

In section 10, we describe how to calculate expectation values of cx-operators. For this 
calculation, one can either transform the expression for the eigenstates of H in the spin representation 
or alternatively, one can transform the expression for the a-operators into the fermionic (Fock) 
representation. We have chosen the second possibility. The transformation from the cx-operators to 
the fermionic operators and is given in eqs. ( f2li| ) and ( |10.3| ) where the (4>^)j are the components 



of the eigenvectors of M defined by eq. ( 2.13; ) (where we use the notation fixed by eq. fl2.9| )). Thus, 



to use this transformation one needs to know the eigenvectors of M. We will now describe how 
to find them analytically in the cases where the zeros of the polynomial are known, following the 
method described in section 3. One first solves eq. ( |3.17| ) to express (pi as function of Tp 1 where the 
coefficients Qij with i,j = 1,2 are given by eqs. (3.18) - (3.21). If x ^ i, the solution for ipi is 
inserted in eqs. ( |3.11| )~( |3.14j ), and the results for the coefficients a,b,f and g used in eq. ( p.8[ ) for 
x 1 respectively eq. ( |3.23| ) for x = 1 yield expressions for ipj and Tpi . The entries of the eigenvector 
(p^ are then given by eq. ( |3.24j ) in terms of ipj and Tpj. In this notation, they still depend on the 
variable x. 

The values of x are obtained as solutions of the polynomial equation (4.2). The polynomial is 



44 



given in the variable z = x 2 . The eigenvectors for the eigenvalues A n with positive real part and 
x n 7^ i are obtained by choosing a square root x n = y/z^ for each zero z n ^ ±1 of the polynomial 
(such that the real part of x n is positive). Observe that due to the quadratic relation between A n 
and x n (eq. (|3.9| )) the values x n and l/x n lead to the same eigenvalue and to the same eigenvector. 

For the eigenvectors corresponding to the eigenvalues A n with negative real part one takes 
x n = —\fzn- The last free parameter Tp l is fixed by the normalization conditions given by eq. 
( 2.20p . Equations ( 2.20Q and ( [2.21 ) are equivalent to the anticommutation relations for the fermionic 



operators. 

For x n = i the equations ( p.3p -(3.7) have to be solved using the ansatz given by eq. ( |3.8| ) for 
(pj and Tpj. Details of this calculation as well as a derivation of the conditions for the appearance of 
additional zero modes on top of the spurious zero mode in the spectrum of Hi ong can be found in 
the appendix. 

13.4- One- and two-point functions of the o~ x -operators 

If no cr z -boundary terms are present in the Hamiltonian or if the condition a__ = a + and /?+ = /?_ is 
met, we obtained formulas for the one-and two-point functions of the aj-operator for both chains 
H and Hi ong . For Hi ong , we considered the ground states given by the eigenstates of ctq and <r£ +1 . 
Remember that Hi ong has a twofold degenerate ground state due to the spurious zero mode. In 
the fermionic language, this corresponds to a vacuum, \vac) and an excited zero mode, |0). Due 
to the symmetry that Hi ong commutes with a x and <r£ , 1 we can pick out the two ground states 



v 



±\ 



vac) ± |0) (see eq. ( |9.3|) ) as eigenstates of a x and cr£ +1 . As far as H is concerned, its ground 
state is given either by \v + ) or by bi owest \v~) where bi owest is the creation operator of the fermion 
corresponding to the energy with the smallest real part as explained above. 

The one-point functions of the crj-operator are non-trivial due to the presence of the non- 
diagonal boundary terms. Without them, they would be zero which is a well-known fact from the 
XX-chain. 

The one- and two-point functions are up to a factor Pfaffians (see eqs. ( |ioT|) , ( [1021) and 



( 10.5| )) of the matrix A given by eq. ( |10.6 ). If no cr z -boundary terms are present, we can further 



reduce this expression to a determinant as given by eq. If a- = a+ and /5+ = /?_, the 

determinant is given by eq. ( |10.14| ). These simplifications are possible because some of the so-called 
basic contractions of pairs of the form given by eq. (10.8) vanish due to the relations (2.40) and (2.42) 



obtained in section 2. To determine the ground state of H, we use the results of this calculation 
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to determine the magnetization (v + \af j+1 \v + ) in section 11 to find the value of the parameter n as 
given by eq. (|11.11| ), 

Our calculation with slight modifications also applies to expectation values of the crj-operator 
with respect to excited states. This is explained in the last paragraph of section 10. 

13.5. List of the results which are going to be used in the following two articles 
Here we give a list of results that we will use in the following two articles: 



Second article [23] 



In order to calculate expectation values of the crj-operator and the <rJ-operator for arbitrary 
position j and lattice length L, we need 

(i) the transformation from the a-operators to the fermionic operators (eqs. ( [2.3; ) and (|10.3j )) 

(ii) the eigenvectors of M (see section 3) and the eigenvalues A n of M (eqs. ([4.1|) and (4.2)) 

(iii) the expressions for the eigenstates of H in the Fock representation (eqs. ( |9,6D , ( |9.7| ) and 
(|9 . 3|) ) and the value of the parameter n defined by eq. (|9.5D 

(iv) the formulas for the one- and two-point correlation functions of crj derived in section 10 

• Third article |4[ 

(i) For the calculation of the excitation spectrum of H in the limit of large L we need the 
polynomial equation (eq. (4.2)) and the projection mechanism (see section 9). 

(ii) For the expressions of the ground state energies in the exactly solvable cases in the limit 
of large L we need the results of table II and III. 

(iii) For the construction of the magnetic charge operator, we need the eigenvectors of M (and 
therewith the eigenvalues, see above). 

14. Discussion 

14-1. Observations on the expressions for the ground state energies 

Up to now we have described how to find the eigenvalues, the ground state energies and the 



eigenvectors for Hi ong (eq. (1.4)) and for H (eq. fll.lD ). We now turn to the discussion of the 
results of our analytical calculations for the cases where the polynomial (eq. (4.2)) can be factorized 
into cyclotomic polynomials. The expressions of the ground state energies of Hi ong and of H are 
given in terms of trigonometric functions only (see tables II and III). It is remarkable that they 
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appear in spite of non-diagonal boundary terms in the Hamiltonians. This reflects the integrability 
of the model. 

Furthermore, notice that for the one-parameter families of exact solutions corresponding to the 
cases 10 - 16 from table II (where the free parameter is called s) the ground state energy has an 
L-independent (but s-dependent) term which appears additively to the L-dependent part of the 
ground state energy. We will show in |2j| that these L-independent contributions to the ground 
state energy of H are related to boundary bound states. 

In some special situations which only appear for non-hermitian Hamiltonians the expressions 
for the ground state energy of H exhibit a rather peculiar behaviour with respect to variations of L. 
Namely, considering one of the cases 10, 14 and 15 and choosing the parameter s in such a way that 
rj = — 1 one observes that for L less than a limiting length Ln m n (which depends on s) the fermionic 
energy 2 Ai owest (which has to be added to the ground state energy of Hi ong to obtain the ground 
state energy of H) is given by twice the -L-independent expression that appears in the ground state 
energy of Hi ong , but with opposite sign (see table III). Therefore this term appears with different sign 
in the expression for the ground state energy of H than in the expression for the ground state energy 
of Hi ong . However, if L is larger than Lu m n, a level crossing in the fermionic spectrum appears and 
another L-dependent fermionic energy becomes smaller than the L-independent energy from before. 
Then this L-dependent fermionic energy has to be added to the ground state of Hi ong instead of 
the L-independent term from before, and the L-independent part no longer switches its sign when 
going from the ground state energy of Hi ong to the one of H . 

Now we discuss the degeneracies in the spectrum of H . Degeneracies may appear due to doubly 
degenerate fermionic energies. In the cases where the polynomial can be factorized into cyclotomic 
polynomials, twofold degenerate fermionic energies can be identified by quadratic factors appearing 
in the factorized form of the polynomial (see table I). - Notice that this observation does not apply 
to the quadratic factors of the form (1 — z) 2 since the polynomial p{z) given by eq. (4.2) has to be 
divided by this term. - In the case where the spectrum of H consists of an odd number of fermionic 
excitations, twofold degenerate fermionic energies also lead to a twofold degenerate ground state. 
The degeneracies in the spectrum of H are also reflected in the partition functions in 24]. 



14-2. Open questions 

Some questions could not be clarified within the framework of this article. 
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a) It is not clear whether table I from section 5 shows all possible cases where the polynomial 
factorizes in six or more factors given by cyclotomic polynomials. Perhaps it is also possible to 
find different factorizations of the polynomial for other boundary parameters which also allow 
to compute all zeros analytically. 

b) From table III one sees that two Hamiltonians may have different boundary terms and still 
have the same spectrum (given by the zeros of the same polynomial) e. g. in the case 2 where 
the Hamiltonian H with a_a + = 1 and (3 + = (3_ = a z = (3 Z = has the same spectrum as 
the Hamiltonian with boundary parameters (3 Z = ^,(*z = 0,a_ = -^e^,a+ = ^e _1 ^,/3+ = 
(3- = 0. This fact gives rise to unknown similarity transformations which remain to be made 
explicit. 

c) For the exactly solvable one-parameter families 10, 14 and 15 with free parameter s we observed 
a surprising behaviour of the expression of the ground state energy of H in the case where the 
parameter rj = — 1. Namely, by increasing the lattice length L and reaching a certain value 
Llimit which is given in terms of s, the L-independent contribution to the ground state energy 
suddenly switches its sign as explained above. The physical origin of this phenomenon is not 
clear. 
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Appendix A. Appearance of fermionic zero modes in the spectrum of Hi ong 

In this appendix, we show how to find the eigenvectors of the matrix M corresponding to the 
eigenvalue zero. This procedure will provide conditions for the boundary parameters which are 
equivalent to the existence of additional zero modes on top of the spurious zero mode. These 
conditions are already mentioned without proof in section 3. 

One might guess that the conditions on the boundary parameters we obtain by constructing 
the eigenvectors are already contained in the polynomial given by eq.(4.2). We will see that this is 
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indeed the case, if we only consider hermitian boundary terms. In the general, non-hermitian case 
this is not true. The polynomial might have more zeros corresponding to eigenvalues A = of the 
matrix M than eigenvectors can be constructed. Therefore, in this case, M is not diagonalizable. 

We are first going to deal with the explicit construction of the eigenvectors. Afterwards we will 
consider the polynomial equation (4.2). 

Appendix A. 1. Construction of eigenvectors 



According to eq.(3.9) A = corresponds to x = ±i. So we solve the boundary equations (3.4)-(3.7) 



using the solution of the bulk equations (|3.3|) given by (|3.8|) with x = i, i.e 



rj "i ; • In i , ? j =g(-i) j + , (A.l) 

where < j < L-\- 1. This can now be used to rewrite the boundary equations in terms of a, b,g, f 
and ipo,Tp 0l (pL + \,Tp L+1 . Introducing the new parameters 

r H75 ±ia ") • r M7! ±i& ) (A ' 2) 

we obtain from the left boundary 

(f = Tp (A.3) 

a_(o-6) = a+(f-g) (A.4) 

r^a + r~b — a + ipo = (A. 5) 

M + »"o/-«-Po = • (A.6) 

The equations from the right boundary give 

(-l) L rpa - r+b + i L+1 /3+^ L+1 = (A.7) 

r+g-(-l) L rpf + i L+1 (3^ L+1 = (A.8) 

(3 + (g+ (-l) L f) = -p_(b+{-l) L a) (A.9) 

fL+l = -Vl+i ■ {A AO) 

Since Tp and tfL+i appear only in the eqs.(A.3) and (A. 10), we have to solve the homogeneous 
system of 6 linear equations given by eqs.(A.4)-(A.9) for the 6 unknowns a, b, f, g, ipo, <fL+x- The 
vector components </? aim <?l+i can then be directly read off from eq.(A.3) respectively eq.(A.lO). 
To have non-trivial solutions for the 6x6 system of equations (A.4)-(A.9) the determinant of 
the corresponding 6x6 matrix has to vanish. This is equivalent to a condition on the boundary 
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parameters aq_, a_, f3 + , i.e. 

a + (3- +a-(3 + = . (A. 11) 

At this point, it is not obvious how many solutions of the eqs.(A.4)-(A.9) may exist. Thus we are 
going to solve them explicitly. In order to do this, we will treat eqs.(A.5)-(A.8) and eqs.(A.4) and 
(A. 9) separately. We first solve eqs.(A.5)-(A.8) for a,b,f,g and then check for consistency with 
eqs.(A.4) and (A. 9). 

Eqs.(A.5) and (A. 7) and eqs.(A.6) and (A. 8) can be rewritten as 

rJ")J t" ),*J 3 )=( L+ r° i. (A.i2) 

where the 2 x 2-matrices R a b and R g f are given by 

R ab =\ rt T ~ a | ,R 9f =\ rt | • (A. 1.1) 

V (-1)% ~r$ J ("1) L+1 ^ 

The determinants of R a b and R g $ have the same values and are given in terms of a z , (5 Z by 



det R a b = det R g f = < 



2/3 z ctz — 1 for L even 

(A.14) 

-V2i(a z + p z ) for L odd 



Once we know the general solution of eq. flA.12j ) , we only have to verify which of them simultaneously 
solve the eqs.(A.4) and (A. 9). Solving eq. (|A.12j ) one has to distinguish two cases: (1) det R a b / , 
(2) det R a b = . Let us first deal with case (1). 



If det R a b 7^ the matrices R a b and R g j can be inverted in order to solve eq.(A.12). Doing this 
we obtain 

a = det ^ - (a+r^tpo ~ i L+1 /3+r~<PL+i) , (A. 15) 

b = ((-l) L a + r^ + i L+1 /?+r> L+1 ) , (A.16) 

9 = ((-l) L a.r^ - i i+1 /3_r^ L+1 ) , (A.17) 

/=deT^(«-> + ii+1 ^a>x +a ) • (A.18) 

Substituting this into eq.(A.4) gives 

(<*+/?_ +a_/3 + )(^ L+1 =0 , (A.19) 
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whereas eq.(A.9) leads to 

(a+/?_ + a-(3 + )ip = 



(A.20) 



Thus, if eq.( A.ll ) is satisfied and if deti? a fe 7^ 0, we obtain two eigenvectors of M corresponding 
to the eigenvalue A = on top of the spurious zero mode because ipo and <pl+i can be chosen 
independently of each other. 

Let us now turn to case (2), i.e. det R a b = 0. Because this condition gives different conditions 
on the non-diagonal boundary terms for L even and L odd respectively (see eq. (|A.14|) ) we discuss 
these cases separately. 

We will first turn to the case where L is odd. Here we have a z = —f3 z which can be read off 
from eq.( A.l4) - Using this in eq.( A.l we can rewrite eq.( A. 12] ) as follows: 

(A.21) 



a + <po = \ L+1 f3+Lp L+ i , a-ipo = -i L+1 f3-ip LA 



a + ip = r f3 a + r^b , a_<^ = r g f + rjg . 



•pi 



For L even we have according to eq,( |A14|) 



a. 



23 



g-. Using this equality, eq.( A.12|) reads 



V2f3 z a + ip = i L f3 + ip L+ i , \/2/3 2 a_ 



■L 



i^/3_^ L+ i 



• z a+ipo 



rJb — r n a 



\V2(5 z a-ifQ = rig -r g f 



(A.22) 



(A.23) 



(A.24) 



Note that at least one of the parameters r„ , rt is different from zero. Thus the eqs. ([A22D and (gjg) 



can be solved either for a and / or for b and g respectively. Note also that due to eqs. ( A.21 ) and 
flOaD ifo and <Pl+i are no longer independent of each other, if one of the parameters a + , a_ , (3 + , /3_ 
is different from zero in contrast to the case (1). 

If all of the parameters a+, a_, /?+, (5- are vanishing the eq. ([A.21 ) respectively eq.( |A.23| ) are 
satisfied automatically. The same holds for eqs. (A. 4) and (A. 9). Thus we simply have to solve 



eqs. (|A.22|) and ( |A.24| ) yielding 



a = (-l) L ^6 , / = (-l) L ^ 5 forr./O 



-\) L \a , g=(-l)Mf forr+/0 

r B r (3 



(A.25) 



(A.26) 



Since on the one hand the parameters a and / or b and g respectively and on the other hand the 
two vector components ifo,<PL+i can be chosen independently we obtain a set of four eigenvectors 
corresponding to the eigenvalue A = on top of the spurious zero mode. 
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If one of the parameters a + , a_ , j3 + , /3_ is different from zero we may solve eqs. ( [A.22 ) and ( [A.24 ) 
for a and / or b and g respectively and use the result in eq.(A.4) to obtain 

aJb = a + g for / (A.27) 

a^a = a + f for r + ^ . (A. 28) 

Additionally using eq.( [A.2l| ) respectively eq.( |A.23 ), we obtain from eq.(A.9) 



= -p_b for / (A.29) 



[3 + f = -0_a for r+ / . (A.30) 



Due to the condition ( |A.11| ), it is always possible to solve eq.( A.27j )-( |A.30| ) and eq.( A.21 ) respectively 



eq.( |A.23 ) by leaving two variables undetermined. The remaining four unknowns can then be given 



in terms of these two. This allows the construction of two further eigenvectors corresponding to the 
eigenvalue A = 0. 

For instance, let us assume a+ ^ 0, ^ and L odd. Eq.( |A.21 ) is then solved by 
ifo = i L+1 ^-ifi + i whereas eqs.( A.27 ) and ( A.2S| ) are solved by g = ^b. The parameters a and / 



a+ 



are then fixed by eq.( A.22| ). Thus a, f,g, ifQ are given in terms of b and <£>l+i which can be chosen 
independently. 

Let us briefly summarize the results of this section. We looked for eigenvectors corresponding to 
the eigenvalue zero. At this point we want to remind the reader that there always exist at least two 
eigenvectors corresponding to the eigenvalue zero, namely the ones which are related to the spurious 
zero mode. We found exactly two further eigenvectors corresponding to the eigenvalue zero of M if 
and only if one of the two following conditions is satisfied : 

(i) a+/?_ + /?+a_ = and j3 z / 7^7 f° r L even respectively a z ^ —f3 z for L odd. 

(ii) a+/3_ + /3+a_ = and at least one of the parameters a±,f3± is different from zero. 
There exist four further eigenvectors corresponding to the eigenvalue zero of M if and only if 
(i) a± = @± = and (3 Z = 7^- for L even respectively a z = —f3 z for L odd. 

There are no other possibilities to have further eigenvectors corresponding to the eigenvalue zero. 
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Appendix A. 2. Zeros of the polynomial at z = — 1 

In this subsection we consider the polynomial equation (4.2) in order to check whether the matrix 
M may have more eigenvalues A = than eigenvectors can be constructed. Since in general M is 
non-hermitian, this may indeed be the case. For hermitian M, i.e. hermitian boundary conditions, 
we will recover the conditions on the boundary terms which we already obtained in the previous 
section. 

Since the polynomial equation is given in terms of the variable z which is related to the 



eigenvalues via eq.(4.1), an additional eigenvalue A = on top of the spurious zero mode corresponds 
to a root of the polynomial at z = — 1. The necessary condition to have at least one eigenvalue A = 
on top of the spurious zero mode therefore translates to 

?(-!) = 

<^D = l + A + B + 2C (A.31) 

^ a-(3 + + a+/3_ = , (A.32) 

where the parameters A, B, C, D, E are defined by eqs.(4.3). These zeros will always appear in pairs 
since with z also - is a zero of q(z). In order to find a condition for the existence of a root at z = — 1 
with higher multiplicity than two, we have to consider the second derivative of q{z). Using eq.(A.31) 
we obtain for even L 

d 2 zP{z)\z=-i = 

2{E - l) 2 + (3 + 2A + B + 2C)L - CL 2 = , (A.33) 

whereas for odd L we get 

l-C-B + AE + (3 + 2A + B + 2C)L - CL 2 = . (A.34) 

Let us now consider eqs.(A.32)-( [A.34 ) for the case of hermitian boundaries. Then eq.(A.32) implies 
that at least one of the parameters a+ = a*_ or /3 + = /?* is equal to zero. Without loss of generality 
we may assume that a+ = 0. This implies immediately C = and 3 + 2 A + B = 2|/3 + | 2 (l + 2a 2 ). 
Using eq.(A.33) we obtain for L even 

2|/3 + | 2 (l + 2« 2 ) = - 2(2a ^- 1)2 . (A.35) 

Since this equality can only be valid if the RHS and the LHS vanish simultaneously, we conclude 
that in the hermitian case further zeros at z = — 1 only exist if 

"+ = /?+ = , a z = — . (A. 36) 
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From eq.( A,34| ) we get for L odd 

W(i + ^ = - l M±^ 

which can only be satisfied if 

a+ = (3+ = , a z = -(5 Z 



(A.37) 



(A.38) 



Further conditions for the existence of more than four zeros at z = —1 can be derived in the same 
manner as eq.(A.33) and eq. (|A.34 ) respectively. Solving eqs.(A.33) and ( A.34|) for A and calculating 
the fourth derivative of p{z) at z = —1 gives the conditions 

4 - 8E + E 2 + 12L(1 - E 2 ) + L 2 (2E 2 - 3B - AC + 8E + 5) + L 4 C = 



for L even and 



C + 8E + B + 6E 2 + 5 + 12L(1 - £ 2 ) + L 2 (7 - 2C - - B + 6£ 2 ) + L 4 C = 



for L odd respectively. Using eqs. QA.36 ) in eq.( A.39| ) gives 

2 



V2p z + 



1 







whereas substitution of eqs.( |A.3§| ) in eq.( A.40| ) yields 



o2_l + L 

z L-l 



(A.39) 



(A.40) 



(A.41) 



(A.42) 



Since neither eq.( [A.41 ) nor eq.( [A.42 ) can be satisfied by any j3 z G R we conclude that in the hermitian 
case we have at most four zeros at z = — 1. It is no surprise that for hermitian boundaries the 
conditions on the boundary parameters obtained in this subsection are equivalent to the ones of the 
previous subsection. However, if M is non-hermitian, the conditions derived in this subsection have 
more solutions than the ones of the previous section. Therefore it may happen that the polynomial 
has more zeros corresponding to an eigenvalue A = than eigenvectors can be constructed. This 
implies that M is non-diagonalizable for certain choices of boundary terms. 
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